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ABSTRACT 

We present a wide-field (4.5 deg 2 ) photometric and spectroscopic survey of the Leo I dwarf spheroidal 
(dSph) galaxy to explore its extended morphology and dynamics. As in previous papers in this series, we take 
advantage of photometry in the M, T2, and DD05 1 filter system to select Leo I red giant branch star candidates, 
and, so far, this selection technique has proven 100% reliable in selecting actual Leo I members among more 
than 100M < 21 .5 Leo I giant candidates having previous or new Keck DEIMOS spectroscopy to a radius > 1 .3 
times the limiting radius of the fitted, central King profile. The two-dimensional distribution of all similarly- 
selected Leo I giant candidates is well fitted by a central single-component King profile of limiting radius 13.3 
arcmin, but many giant stars are found outside this newly derived King limiting radius. The density profile thus 
shows a break at a major axis radial distance of ~ 10 arcmin produced by an excess of stars at and beyond the 
King limiting radius (spectroscopically confirmed to be made of true Leo I members), and primarily along the 
major axis of the main body of the rather elongated satellite. This spatial configuration, a rather flat velocity 
dispersion profile and an asymmetric radial velocity (RV) distribution among the Leo I members at large radii 
together support a picture where Leo I has been tidally disrupted on at least one, but at most two, perigalactic 
passages of a massive Local Group member. We demonstrate this hypothesis using mass-follows-light, /V-body 
simulations of satellites in a Milky Way-like potential that reproduce the observed structural and dynamical 
properties of Leo I remarkably well. These models include ~3x 10 7 solar mass, tidally disrupting dSph 
systems on bound orbits with rather high eccentricity (0.93-0.96) and small perigalactica (10-15 kpc). The 
simulations allow the first observationally constrained orbit for Leo I without the measurement of its proper 
motion and show that the observed RV distribution is more consistent with a two Milky Way orbit history for 
the satellite while ruling out a Leo I orbit that includes a previous association with M31 within the last 10 Gyr. 
Given the overall success of tidally disrupting mass-follows-light satellite models to account for the observed 
properties of Leo I, we conclude that there is no need to invoke an extended dark matter halo around the satellite 
(e.g., as one explanation of the velocity dispersion and radial profiles at large radii), and that an overall modest 
M/L for the satellite is consistent with the available data. That a satellite on such a large (apogalacticon of 
~ 450 kpc), long period (P ~ 6 Gyr) orbit as Leo I can experience tidal disruption suggests that similarly 
structured satellites with even smaller (eccentric) orbits will experience even greater tidally-induced mass loss 
rates. 

Subject headings: galaxies:evolution — galaxies:interactions — galaxies:halos — galaxies:individual(Leo I) 
— galaxies:photometry — galaxies:structure 
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1. introduction 

1.1. Motivations for a New Study of the Leo I System 

Modern theories of the evolution of structure in the Uni- 
verse that include cold dark matter (CDM) grow galaxies and 
clusters of galaxies and their dark halos through the accumu- 
lation of smaller subunits (e.g., White & Rees 1978; Davis et 
al. 1985; Navarro, Frenk, & White 1996, 1997; Moore et al. 
1999). But the notion that the Milky Way (MW) halo was 
built up by protracted infall of "protogalactic fragments" af- 
ter an initial central collapse had already been established by 
Searle & Zinn (1978) based purely on stellar population ar- 
guments. In CDM scenarios dwarf satellite galaxies represent 
the visible parts of (predominantly dark) subhalos; however, 
the number of dwarf galaxies discovered so far is several or- 
ders of magnitude less than the predicted number of subhalos 
made by CDM simulations (Kauffman, White, & Guiderdoni 
1993; Klypin et al. 1999; Moore et al. 1999). This may be 
an indication that the majority of (especially smaller) subha- 
los have either not formed stars (Bullock, Kravstov, & Wein- 
berg 2001) or have been destroyed (Hayashi et al. 2003), and 
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that the visible satellites of today represent the high mass end 
of the mass spectrum of DM subhalos (Stoehr et al. 2002; 
Hayashi et al. 2003) or those that were able to accrete substan- 
tial amounts of gas before reionization (Bullock, Kravtsov, & 
Weinberg 2000; Benson et al. 2002; Taylor et al. 2004). 

On the other hand, dwarf galaxies exhibit some proper- 
ties that may be inconsistent with the expected properties of 
CDM subhalos. For example, the flat central density pro- 
files of dwarf spheroidal (dSph) galaxies are at odds with the 
cuspy interiors predicted by CDM (e.g., Navarro, Frenk, & 
White 1996, 1997; Moore et al. 1998), although this might 
be reconciled by appealing to triaxial halos (Lokas 2002; 
Navarro 2004). Alternatively, warm dark matter (WDM) al- 
lows smaller central phase-space densities, and studies of 
dSph systems have provided important constraints on the 
properties of WDM species (Lin & Faber 1983; Gerhard & 
Spergel 1992; Goerdt et al. 2006; Strigari et al. 2006). In 
addition to the problem of the central concentrations of satel- 
lites, their apparent alignments around parent halos — as has 
been observed in the Milky Way and argued to be related to 
"dynamical families" of satellites (Kunkel 1979; Lynden-Bell 
1982; Majewski 1994; Fusi Pecci et al. 1995; Lynden-Bell 
& Lynden-Bell 1995; Palma, Majewski, & Johnston 2002) — 
have also been used to question the viability of dwarf galaxies 
as putative subhalos (Kroupa, Theis, & Boily 2005; Kang et 
al. 2005). On the other hand, such alignments have also been 
reported in high resolution N-body, hierarchical structure sim- 
ulations, where subhalos accrete along filaments (Libeskind et 
al. 2005; Zentner et al. 2005; Wang et al. 2005; Libeskind et 
al. 2007). 

Despite uncertainties over the precise connection of ob- 
served satellite galaxies to the prevailing theoretical CDM 
models, it has become observationally clear that satellite 
galaxies (e.g., Ibata, Gilmore, & Irwin 1995; Newberg et 
al. 2002; Majewski et al. 2003; Yanny et al. 2003; Crane 
et al. 2003; Rocha-Pinto et al. 2003) and even star clusters 
(Odenkirchen et al. 2001, 2003; Rockosi et al. 2002; Lee et 
al. 2004; Grillmair & Johnson 2006; Belokurov et al. 2006c; 
Grillmair & Dionatos 2006) could be significant contribu- 
tors to the luminous halo of the MW and other galaxies. N- 
body simulations (Oh, Lin, & Aarseth 1995; Piatek & Pry or 
1995; Johnston, Spergel, & Hernquist 1995; Johnston, Sig- 
urdsson, & Hernquist 1999; Johnston, Choi, & Guhathakurta 
2002; Munoz, Majewski, & Johnston 2007) show how a dwarf 
galaxy can experience tidal disruption in its outer parts during 
close encounters with the central potential well of MW-like, 
parent systems. Stars that escape the system will form tidal 
tails, like those observed in the Sagittarius (Sgr) system (e.g., 
Ibata, Gilmore, & Irwin 1995; Newberg et al. 2002; Majew- 
ski et al. 2003; Belokurov et al. 2006d). Such extratidal stars 
can create "breaks" in the projected radial star-count profiles 
of the satellites as the density law transitions from a steeply- 
declining central density law for bound stars to a much more 
gradual decline at radii where unbound stars start to contribute 
significantly. Such "breaks" have been observed in several 
MW dwarf spheroidal (dSph) satellites (Eskridge 1988a,b; 
Irwin & Hatzidimitriou 1995; Majewski et al. 2000a, 2003; 
Palma et al. 2003; Wilkinson et al. 2004; Westfall et al. 2006; 
Siegel et al. 2005; Munoz et al. 2006). For example, the ra- 
dial density profiles of nearly every one of the then-known 
dSphs of the Milky Way in the thorough Irwin & Hatzidim- 
itriou (1995, hereafter IH95) study (see their Figure 2) show 
an excess of stars with respect to the outer parts of their best- 
fitting model King profiles. In the particular case of Leo I, 



because of its extreme distance, the number of dSph stars rel- 
ative to the background level in the outermost regions was 
too low to judge conclusively the significance of its apparent 
break and limited the ability of IH95 to explore the Leo I pro- 
file significantly past the King limiting radius with their pho- 
tographic star counts. Undoubtedly, a radial surface density 
profile with better signal-to-noise and spatial coverage will 
help verify whether this extremely distant Galactic compan- 
ion shows the break profile trait seen in other MW dSphs. 
Determining the structure of Leo I to large radii is a primary 
goal of the present study. 

A second goal is to increase our understanding of the inter- 
nal dynamics of Leo I. If dSph galaxies are the visible coun- 
terparts of the largest DM lumps (e.g. Stoehr et al. 2002), one 
might expect the process of tidal disruption to be inhibited 
or lessened in these massive subhalos. On the other hand, if 
dSphs can be proven to have tidal tails, the latter would place 
constraints on the dark matter components within these ob- 
jects (Moore 1996). The Sgr dSph provides a vivid example 
of a tidally disrupting system for which the degree of disrup- 
tion, as measured by the nature of the tails, has been used to 
constrain its dark matter content (Law, Johnston, & Majew- 
ski 2005); on the other hand, the Sgr system is sometimes 
considered to be an exception to the norm for dSph galaxies 
(Mateo et al. 1998; Mateo 1998). That Sgr may not be an ex- 
ceptional case is suggested by studies of possible tidal disrup- 
tion in other MW dSphs (Gould et al. 1992; Kuhn, Smith, & 
Hawley 1996; Majewski et al. 2000a, 2006; Gomez-Flechoso 
& Martmez-Delgado 2003; Palma et al. 2003; Walcher et al. 
2003; Westfall et al. 2006; Munoz et al. 2005; Munoz et al. 
2006; Munoz, Majewski, & Johnston 2007). 

The effects of tidal disruption can be inferred not only by 
the existence of "break" in the density profiles of these other 
MW dSphs, but also, more recently, by the velocity charac- 
teristics of the dSphs at large radii. The past decade has seen 
substantial progress in measuring radial velocities for large 
numbers of stars in at least the more nearby MW dSph galax- 
ies (Armandroff et al. 1995; Tolstoy et al. 2004; Wilkinson 
et al. 2004; Munoz et al. 2005; Munoz et al. 2006; Walker et 
al. 2006). These studies not only imply large central mass-to- 
light ratios in many of the satellites, but generally reveal rela- 
tively flat velocity dispersion profiles reaching into the break 
profile regions (Munoz et al. 2005; Munoz et al. 2006), a phe- 
nomenon reproducible in N-body simulations of disrupting, 
mass-follows-light satellites (Munoz, Majewski, & Johnston 
2007). 

However, a tidal disruption interpretation of these velocity 
dispersion profiles to large radii is not unique based on the ex- 
tant data, because flat velocity dispersion profiles may also be 
accommodated to arbitrary radii by extending the DM halos 
in which the luminous dSphs are embedded. These "equilib- 
rium models" imply mass-to-light ratios that increase with ra- 
dius and boost the global ratio — in some cases quite substan- 
tially (e.g., Kleyna et al. 2002). Knowing whether mass fol- 
lows light in dSphs or if their luminous components lie within 
extended dark halos is critical to establishing the regulatory 
mechanisms that have inhibited the formation of galaxies in 
all subhalos and the extent to which luminous satellites are 
vulnerable to disruption and populate the stellar halos of L* 
galaxies. 

Leo I provides an interesting contrasting case to other 
dSphs in that (1) this particular satellite clearly has an un- 
usual orbit that has protected it from but a few potential tidal 
encounters, and (2) unlike the apparently high M/L systems 
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(like Ursa Minor, Draco and Carina), the previously measured 
mass-to-light ratio of Leo I is only ~ 6Mo/Loy, i.e. at the 
low end of the dSph M/L scale. Of course, it is of interest to 
know whether and how these two traits of this satellite may re- 
late, and whether the properties of Leo I, if explored more ex- 
tensively, may lend new insights into the issues raised above. 
For example, as we shall show in this paper, new signatures 
found in the velocity distribution of Leo I stars at large radii 
(namely, an asymmetry) may provide a way to break the de- 
generacy in interpretations of the velocity dispersion profile 
for at least some dSphs. 

The heart of our study is a new photometric survey cover- 
ing 4.5 deg 2 of the sky centered on Leo I, and new Keck spec- 
troscopy of 105 more Leo I giant stars than available in the 
literature. As with other contributions in this series, we adopt 
a technique based on multicolor Washington filter imaging, 
including the use of the DD05 1 filter centered on the gravity- 
sensitive MgH+Mgb triplet spectral feature near 5150 A to 
identify giant stars associated with the Leo I dSph. We use 
bright K giants of Leo I as a tracer population of the struc- 
ture of the dSph because (1) they are the most easily detected 
type of stars over large areas with the use of a mid-size tele- 
scope, (2) the Washington+DD(95 1 separation technique has 
already been proven successful in the discrimination of metal- 
poor dSph giant stars from foreground metal-rich field dwarf 
stars and thereby to significantly increase the signal-to-noise 
of the faint, diffuse outer parts of Local Group dSphs (e.g., 
Majewski et al. 2000a; Palma et al. 2003; Munoz et al. 2005; 
Munoz et al. 2006), and (3) these stars are accessible to spec- 
troscopic membership and dynamical follow-up with the cur- 
rently largest telescopes. 

Because we are able to explore the Leo I dSph to large radii 
with much better signal-to-noise than previous studies, we can 
rederive the structural parameters for Leo I (§4.1), which are 
important to assessing the M/L using the traditional, King 
(1962) methodology based on the internal dynamics of the 
system. A newly derived M/L for Leo I in this way is dis- 
cussed in §4.2. We also explore the possible mass distribu- 
tion for Leo I within the context of a tidal disruption scenario 
using iV-body simulations of Leo I-like satellites orbiting a 
MW-like galaxy in §6.2. 

The unusually high radial velocity of Leo I at its extreme 
distance naturally leads to interesting questions about its spe- 
cific orbit. There have been two studies that are in disagree- 
ment about potential orbits for Leo I. On one hand, Byrd et al. 
(1994) conclude that Leo I was once loosely bound to M31 
and now is in an unbound, hyperbolic orbit about the MW. 
They estimate Leo I's MW perigalactic distance to be <~ 70 
kpc and for this nearest approach to have occurred 2-4 Gyr 
ago. On the other hand, in an earlier study, Z89 concluded 
that Leo I probably did not originate in the M3 1 system and 
that the most reasonable assumption is that Leo I is bound to 
the MW. An unexpected benefit of our study is that we have 
been able to derive new constraints on the Leo I orbit from 
detailed study of its structure and dynamics (§6.2), and we 
find reasonable agreement of these results to inferences about 
Leo I's orbit based on its star formation history. 

Dwarf spheroidal galaxies, with other halo objects such as 
globular clusters and field stars, are also useful test particles 
for probing the large scale mass distribution and total mass 
of our Galaxy. Despite much work, the spatial extent and to- 
tal mass of the MW remain among the more poorly estab- 
lished Galactic parameters. While traditionally thought to be 



a "lesser sibling" to the Andromeda galaxy (M31) in terms of 
mass, some recent work (Cote et al. 2000; Evans & Wilkin- 
son 2000; Evans et al. 2000; Geehan et al. 2006; Seigar et 
al. 2007) suggests that, in fact, the MW edges out M31 as 
the most massive galaxy in the Local Group. With the tradi- 
tional analyses, the mass estimate is sensitive to the inclusion 
or exclusion of one particular object, Leo I because of its com- 
bined unusually large Galactocentric radial velocity (+177 ± 3 
km/s; Zaritsky et al. 1989, hereafter Z89) and great distance 
(257 ±8 kpc, see §1.2). Estimates of the Galactic mass by 
Z89 varied by a factor of 3-4 depending on whether Leo I is 
considered as bound to the MW or not. While recent stud- 
ies by Wilkinson & Evans (1999), Sakamoto, Chiba, & Beers 
(2003) have decreased the dependency of MW mass estimates 
on including Leo I by using Bayesian likelihood methods and 
larger samples of halo objects, Leo I is still considered to be a 
determining factor for fine tuning the results. 

1.2. Previous Photometric Studies of Leo I 

The Leo I dSph was discovered more than half a century 
ago by Harrington & Wilson (1950) during the first Palomar 
sky survey. Due to its extreme distance and its angular prox- 
imity to the 1 st-magnitude foreground star Regulus, photo- 
metric studies of Leo I have been difficult. Not until CCD ar- 
rays were developed were the first color-magnitude diagrams 
(CMD) of Leo I constructed. Fox & Pritchet (1987), Reid 
& Mould (1991), Lee et al. (1993), and Demers, Irwin, & 
Gambu (1994) presented early ground-based CMDs. Later, 
Caputo et al. (1999) and Gallart et al. (1999a) presented Leo I 
CMDs based on data taken with the Hubble Space Telescope 
(HST) and they reach the main- sequence turnoff (MSTO) of 
the oldest (> 10 Gyr) Leo I populations. These deep CMDs 
allowed detailed studies of the multiple stellar populations 
and the complex star formation history of Leo I (Gallart et 
al. 1999b). Leo I was thought to be unique among the MW 
satellite dSphs for not having a conspicuous horizontal branch 
(HB) population until a 12'x 12' ground-based survey on 
Leo I by Held et al. (2000) revealed an extended HB struc- 
ture in its CMD. More recently, Held et al. (2001) discovered 
more than 70 RR Lyrae variables with pulsational properties 
suggesting an intermediate Oosterhoff type similar to other 
dwarf galaxies in the Local Group (Siegel & Majewski 2000; 
Cseresnjes 2001; Priztl et al. 2002). The existence of both ex- 
tended blue HB population and RR Lyrae stars suggest that 
Leo I is in fact similar to other local dSph galaxies in having 
a > 10 Gyr population, likely formed in the initial collapse 
of the system. The extended star formation history of Leo I, 
which includes this initial starburst, followed by a quiescent 
phase and a new burst of star formation activity starting ~ 7 
Gyr ago (Gallart et al. 1999b) may be intimately tied to its or- 
bital dynamics, since close interactions between galaxies are 
known to be triggers of star formation. 

The distance to Leo I has been derived using various pho- 
tometric methods. Lee et al. (1993) used the tip of red giant 
branch (TRGB) method to derive the Leo I distance modu- 
lus (m-M)o = 22.18 ± 0.1 1, while Demers, Irwin, & Gambu 
(1994) used both the apparent magnitudes of red clump and 
of the carbon stars to estimate (m-M)o = 21.56 ±0.25. Held 
et al. (2001) used the mean magnitude of RR Lyrae variables 
in Leo I and derived (m-M)o = 22.04 ±0.14, and more re- 
cently, Bellazzini et al. (2004) provided a new estimate of 
(m-M)o = 22.02 ±0.1 3 using the TRGB method. A weighted 
average of the distance modulus derived in the four studies 
gives (m-M)o = 22.05 ±0.07, which converts to a distance of 
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257 ± 8 kpc. Throughout this study, we adopt these values. 

The following sections include a presentation of the data 
from our photometric survey of Leo I (§2), and a description 
of the photometric identification of the Leo I giant star can- 
didates (§3). The two-dimensional distribution of Leo I gi- 
ant star candidates is discussed in §4, and new morphological 
parameters for the dSph are derived. In §5 we present new 
Keck Observatory spectroscopy of a subsample of our Leo I 
giant candidates and in §6 we discuss the implications of our 
results, making use of new A^-body simulations of a tidally 
disrupting dSph satellite galaxy that appear to generate sim- 
ilar Leo I properties to those we have observed. Finally, a 
summary of our work and conclusions are in §7. 

2. PHOTOMETRIC OBSERVATIONS AND DATA REDUCTION 

The images used in this study were obtained with the 
Mayall 4 meter Telescope at Kitt Peak National Observa- 
tory (KPNO) during the nights of UT 1998 November 17 
and 2002 May 2-5. We used the Mosaic I 8Kx8K CCD, 
which has a pixel scale of 0"26 per pixel resulting in a 
36 x 36 arcmin 2 field of view. The camera is an array of eight 
2048 x 4096 CCD chips. We used the broad-band Washington 
M, Harris /, and intermediate-band DD051 (hereafter, D51) 
filters. 11 Our survey fields were selected to lie predominantly 
along the major axis of Leo I. The area to the north of Leo I 
was sampled but a large part to the south of Leo I was avoided 
due to Regulus. All program fields were overlapped by 6' with 
adjacent fields as a check on consistency of the photometry. 
Table 1 summarizes the basic information of the CCD fields 
used in this study. 

The raw images were pre-processed using the CCDPROC 
task in the IRAF MSCRED package. 12 The flat-fielding was 
done with special care since the KPNO T2 and D5 1 passband 
images are affected by a pupil image that produces an artifact 
that increases the amount of background light near the center. 
The T2 and D5 1 images also suffer from fringing. We care- 
fully followed the procedures described in Valdes (1998) to 
correct for these effects. Once the pre-processing was done, 
we split each Mosaic image into its eight sub-images and per- 
formed stellar photometry via the DAOPHOT II/ALLSTAR 
(Stetson 1987) package. A point-spread function (PSF) was 
constructed using 25-50 bright and isolated stars for each sub- 
image. The quality of each PSF was improved by removing 
neighboring stars and reconstructing the PSF iteratively. PSF 
magnitudes were derived using UNIX shell scripts based on 
ALLSTAR. The growth-curve analysis package DAOGROW 
(Stetson 1990) was then used to correct for the missing light 
lying outside of the PSF tail (the aperture correction). 

Measured instrumental magnitudes were calibrated against 
Geisler ( 1 990, 1 996) standards that were observed many times 
in different airmass ranges over each observing run. We fit a 
transformation equation of the form: 

MAG-mag = a\ + a2X + a^C, (1) 

where MAG is the Geisler standard magnitude, mag is the in- 
strumental magnitude, X is the observed airmass, and C is 
the color index. We tested for terms in XC and these were 
found to be negligible. The colors M-T2 and M-D51 were 

1 ' Hereafter, we denote the / filter as T2 filter since their response curves 
are nearly identical (see, e.g., Fig. 1 of Lejeune & Buser 1996 and discussion 
by Majewski et al. 2000a). 

12 IRAF (the Image Reduction and Analysis Facility) is distributed by the 
National Optical Astronomy Observatories. 



used in the equations derived for the M and D5 1 filters, re- 
spectively. A color term for the T2 filter was also found to be 
unnecessary. The RMS of the solutions for the transformation 
equations were less than 0.01 magnitude in all three filters. 
For the 1998 observations, we used the existing transforma- 
tions derived for other observations on these same nights by 
Ostheimer (2002). He notes that on the night of UT 1998 
November 17 there were minor transparency variations, so 
this night's data are considered non-photometric. We subse- 
quently tied the 1998 photometry to that of May 2002 in the 
following manner: First, multiply-measured stars in the over- 
lapping regions among the 2002 observations were used to 
derive and apply frame to frame offsets, which in most cases 
were less than 0.01 magnitudes for all three filters. This en- 
sures that all of the 2002 observations share the same photo- 
metric zero-point. We then used stars in overlapping regions 
of the 1998 and 2002 observations to calculate the average 
magnitude offsets of each non-photometric frame relative to 
the photometric ones, and the corresponding offsets were ap- 
plied to all objects in each non-photometric frame. We it- 
erated these steps until all of the average offsets among all 
frames were less than 0.001 magnitudes. 

We note that images taken with mosaic CCDs such as those 
used in this study may suffer from chip-to-chip sensitivity dif- 
ferences that could make the color terms be chip-dependent. 
Ostheimer (2002) found that the chip-to-chip color terms of 
the KPNO Mosaic to be negligible by observing the standard 
fields in Washington and DD05 1 filters on every CCD chip in 
the Mosaic array and by cross-comparing the standard coeffi- 
cients for each chip. 

Astro me try of the detected objects was obtained by run- 
ning the TFINDER task in the IRAF FINDER package and 
using USNO-A V2.0 catalog stars (Monet et al. 1998) as ref- 
erence. Since Leo I is located at a high Galactic latitude 
(b = 49°), the foreground reddening is not significant [accord- 
ing to Schlegel, Finkbeiner, & Davis 1998, E(B-V) = 0.037]. 
It is, however, important to consider the variation of reddening 
within and across our program fields which may result in sys- 
tematic differences in stellar magnitudes and colors. Each ob- 
ject in our data set has been corrected for reddening based on 
its Galactic coordinates (converted from the equatorial coor- 
dinates) and direct reference to the reddening map constructed 
by Schlegel, Finkbeiner, & Davis (1998). The E(B-V) of all 
objects in our fields range from 0.03 1 to 0.050. 

Our survey covers a large spatial area and our photometry 
is likely to contain a large number of galaxies. Bad columns, 
random cosmic rays, and photoelectron bleeding from satu- 
rated stars are other possible nonstellar contaminants. Non- 
stellar objects were eliminated from our photometry using the 
two DAOPHOT image quality diagnostic parameters SHARP 
and x- All sources beyond the range -0.3 < SHARP < 0.3 
were considered as objects with extreme non-stellar morphol- 
ogy, and were rejected. In case of x 13 , our experiments show 
that the acceptable range changes as a function of radial dis- 
tance from the center of Leo I because crowding near the cen- 
ter affects the stellar profiles. Therefore, we have elected to 
use a different \ selection criteria for sources inside a major 
axis radial distance 14 of 8 arcmin (x< 2.5) and outside this 

13 The ratio of the observed pixel-to-pixel scatter from the image profile 
compared to the expected pixel-to-pixel scatter from a model stellar image 
profile. 

14 Until we derive new structural parameters of Leo I using our photomet- 
ric data, we use those of IH95, i.e. r t = 12'.6, PA = 79°, e — 0.21, for defining 
the tidal boundary of Leo I. 
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radius (x< 1.3). 

In Figure 1, we plot the DAOPHOT internal photometric 
error for the final catalog stars as functions of calibrated mag- 
nitudes in the three different bands. 

Figure 2a shows the sky distribution of all stars detected by 
DAOPHOT in celestial coordinates. The varying density of 
sources reflects differences in limiting magnitude and seeing 
among our different fields. Figure 2b shows the distribution 
of sources to a uniform depth. In Figure 3, we show the (M- 
T 2 ,M) and (M- T 2 , T 2 ) Q CMDs for the stars in Figure 2a. The 
left and right panels of Figure 3 show the CMDs for stars 
in the central field of Leo I (C field) and in all other fields, 
respectively. The dominant CMD structure seen in the left 
panels is the upper part of Leo I red giant branch (RGB). The 
prominent Leo I red giant clump (RGC) is present at 21.5 < 
M < 23.0 (20.0 < T 2fi < 22.0) and 0.6 < (M-T 2 ) < 1.6. A 
small clump of asymptotic giant branch (AGB) stars at Mo ~ 
19.5 (T 2fl ~ 17.5) and (M-T 2 ) Q ~ 2.0 is also apparent. One 
other noticeable feature in the left panels of Figure 3 is the 
group of stars at 0.0 < (AT- T 2 ) Q < 0.7 and 19.0 < M < 21 .5 
(19.5 < T 2fi < 21 .0). These stars are thought to be anomalous 
Cepheids or short period Cepheids of a few hundred Myr age 
(Gallart et al. 1999a). The right panels of Figure 3 are more or 
less typical CMDs for a high Galactic latitude field containing 
a mixture of Galactic disk stars, halo giants, field HB and blue 
straggler stars. The sharp edge at (M- T 2 )o ~ 0.65 is due to 
the main-sequence turnoff (MSTO) of field stars smeared in 
apparent magnitude by the range of distances along the line 
of sight (see e.g., Reid & Majewski 1993; Chen et al. 2001; 
Siegel et al. 2002). The many stars bluer than (M- T 2 ) Q = 0.65 
are likely field HB stars and blue stragglers. However, some 
low density of Leo I stars of all types may also lie in the outer 
fields, swamped by the Milky Way foreground. 

3. IDENTIFICATION OF LEO I GIANT STAR CANDIDATES 

The methodology we use in this study to select clean sam- 
ples of Leo I giant star candidates is adopted from that used 
by Majewski et al. (2000b) and subsequent papers in this se- 
ries. In summary, we photometrically select stars that have 
(1) magnesium (MgH + Mgb) band/line strengths consistent 
with those of giant stars, and (2) combinations of surface tem- 
perature and apparent magnitude consistent with the RGB of 
Leo I. In the following sections, we demonstrate the applica- 
tions of these two criteria using our M, T 2 , and D51 photom- 
etry. To improve the reliability of our Leo I giant star catalog, 
we restrict our sample to the stars that have photometric errors 
less than 0.04 mag in each band. 

3.1. Giant Star Discrimination in the Color-Color Diagram 

The D51 filter measures the strength of the MgH+MgZ? 
spectral feature which is a good indicator of stellar surface 
gravity. Since dwarf and giant stars are differentiated by their 
surface gravities, we use the (M-D51) colors — where M 
filter acts as a continuum measure against D5 1 — combined 
with the surface temperature sensitive (M-T 2 ) color 15 to sep- 
arate the foreground dwarf stars and distant metal-poor giant 
stars. The use of the (M-T 2 ,M-D51) color-color diagram 
as a dwarf/giant separating tool has been utilized in several 
studies (Majewski et al. 2000a,b; Morrison et al. 2000; Palma 
et al. 2003; Westfall et al. 2006; Munoz et al. 2005; Munoz et 
al. 2006). 

15 Majewski et al. (2000a) have shown that the (M—T2) has a linear re- 
lationship with V -I , which is a good surface temperature indicator for late 
type stars. 



Figure 4 shows the color-color diagrams for all stars that 
survived the magnitude error, SHARP, and \ rejection. To 
better show the difference of the (M -T 2 ,M - D51) distribu- 
tion between two distinct components, i.e. Leo I giants and 
foreground dwarf stars, we divide our plot into stars in (a) the 
C field alone and (b) all other fields. The central part of the top 
panel of Figure 4 is dominated by Leo I giant stars whereas the 
lower panel shows mainly the elbow-shaped, high-metallicity 
dwarf star locus and only a few stars in the classic "giant" re- 
gion of the diagram (see Majewski et al. 2000b). The solid 
bounding box is defined as follows. The bottom and right 
boundaries were set by using the distribution of Leo I giant 
stars in the top panel. The lower-left boundary was drawn 
roughly parallel to the dwarf locus but offset by about +0. 1 
magnitude in M - D5 1 to account for the increased color er- 
rors at the faint end of our data set. We note from Figure 4a 
that this conservative limit does result in the loss of the bluest 
Leo I giant stars, but ensures that few dwarf stars will scat- 
ter into our selection. We truncate the upper-left boundary 
of the box at (M - T 2 )o =1.0 since this corresponds to the 
blue boundary of our color-magnitude selection region dis- 
cussed below. Stars are considered to be giant candidates if 
they are inside the bounding limits in the two-color diagram. 
Our color-color selection rejects all but the most metal-poor 
([Fe/H] < -2.0) foreground dwarf stars, whereas most dwarf 
stars at higher metallicities lie mainly on the elbow-shaped 
locus (Paltoglou & Bell 1994; Majewski et al. 2000a). 

3.2. Identifying Leo I Giant Star Candidates in the 

Color-Magnitude Diagram 

To select a purer sample of Leo I giants, we now apply a 
second, color-magnitude criterion to our giant candidate list. 
Regardless of the angular separation from the core of Leo I, 
any giant star in Figure 2 associated with Leo I should have 
a combination of (M- T 2 ) and M magnitude that places them 
within the giant branch of Leo I. The left panel of Figure 5 
shows the CMD for giant candidate stars within the IH95 tidal 
boundary of Leo I, while the right panel is for all of the giant 
candidate stars. We use the stars in the left panel to delineate 
the RGB structure in the CMD. The boundaries were drawn 
to roughly follow the overall structure of the giant branch. 
The box is slightly extended to a point brighter and redder 
than the tip of RGB to include some Leo I stars on the tip of 
the asymptotic giant branch (AGB). The lower boundary was 
set at Mo = 21.5 where the luminosity function of our color- 
color selected stars turns overdue to incompleteness. We note 
that the nine stars brighter than Mo = 21 .5 lying outside of the 
RGB selection region in Figure 5a are likely other (asymp- 
totic) giant members of Leo I, but we elect to exclude these 
few stars in the interests of maintaining a conservative selec- 
tion that excludes interloping non-members. This CMD se- 
lection criterion is finally applied to the entire sample of giant 
candidates as shown in Figure 5b. In this way, a total of 1282 
stars are selected as Leo I giant candidates. 

3.3. Leo I Giant Stars in Earlier Spectroscopic and 

Photometric Studies 

It is worthwhile to compare the results of our giant star 
selection criteria against known spectroscopically and photo- 
metrically confirmed Leo I members. M98 obtained spectra 
of 32 stars within the core radius (r c ) of Leo I and another one 
just outside of r c . The results from their derived RVs show 
that all 33 of these stars in this densest part of Leo I are mem- 
bers of the dSph. There have also been carbon star surveys 
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of Leo I using spectroscopic (Azzopardi, Lequeux, & Wester- 
lund 1985, 1986) and photometric data (Demers & Battinelli 
2002). Stars 21 and 24 of M98 sample are the same stars listed 
as carbon stars Azzopardi 18 and Azzopardi 1, respectively in 
Azzopardi, Lequeux, & Westerlund (1985, 1986). We have 
cross- identified 28 red giant stars of M98, 14 carbon stars of 
Azzopardi, Lequeux, & Westerlund (1985, 1986), and 2 car- 
bon stars of Demers & Battinelli (2002) with our photometry. 
Other stars in the past literature that were not cross-identified 
in our survey are missing because they are situated within the 
physical gaps of the Mosaic CCD chips. 

Figure 6 shows the identified giant (filled squares) and car- 
bon stars (filled triangles) in the color- color diagram (top 
panel) and the CMD (lower panel). The two reddest carbon 
stars (inverted filled triangles) are from Demers & Battinelli 
(2002). It is immediately apparent that the carbon stars are 
not situated within the giant star selection region in the color- 
color diagram. Instead, these stars extend to the lower right 
part of the diagram indicating that the carbon stars have strong 
absorption features within the D5 1 filter bandwidth. In fact, 
carbon stars are well known to have a strong C2 Swan band 
absorption feature at a bandhead of A 5 165 A (see e.g., Fig. 
1 of Christlieb et al. 2001, for a representative spectrum of 
an R-type carbon star), well within the bandwidth of the D5 1 
filter. Despite the weak surface gravity of these cool giant 
stars, the strong C2 absorption lines dominate the M-D51 
color and place the carbon stars away from the giant selection 
regions shown in Figure 6. On the other hand, all but three 
of the normal giant stars in the M98 sample are within the 
color- color selection box (Figure 6a). We could include these 
three stars by extending the lower boundary of our color-color 
selection box, but because it might introduce more contami- 
nants (such as metal-poor subdwarfs from the thick disk or 
halo), we opted not to do so. Our RGB selection in the lower 
panel includes all M98 giant stars. The placement of filled 
squares in Figure 6 thus substantiates the reliability of our 
selection technique, though we note that all cross-identified 
non-carbon giant stars of Leo I lie in the bright end of the gi- 
ant branch where photometric errors are on average only 0.01 
magnitude in all three bands. However, in §5 we present new 
spectroscopy showing that the catalog maintains its reliability 
to significantly fainter magnitudes. 

3.4. Spatial Distribution of Leo I Giant Star Candidates 

Figure 7 shows the sky distribution of all stars selected as 
Leo I RGB stars by our color-color and color-magnitude se- 
lection methods. Of particular interest are the 52 stars that lie 
outside the limiting radius found by IH95 and 26 stars out- 
side twice the IH95 limiting radius. The inner part of the 
spatial distribution (lower panel) shows more clearly the fall 
off in the concentration of Leo I giant candidates. Interest- 
ingly, we find that stars spill over the IH95 boundary from the 
main body of Leo I predominantly along the semi-major axis 
to both sides. We discuss the more extended distribution of 
"Leo I giant candidates" below. 

3.5. Evaluation of Background Level 

The two-criteria Leo I giant star selection method described 
above is designed to maximize the reliability of our candidate 
selection. However, our final sample may still contain con- 
taminants. Any halo giant or extreme ([Fe/H]< -2.0) sub- 
dwarf with an observed (M-T2,M)o combination similar to 
our Leo I RGB stars can be caught in our color-color and 



color-magnitude selection net. So too can normal metallic - 
ity dwarfs artificially scattered into the selection from photo- 
metric errors (these will be main sequence turnoff dwarfs or 
subgiants with My < 6, so that these will also be predomi- 
nantly from the Milky Way halo). The number of these con- 
taminants can be estimated by a simple technique described 
in Majewski et al. (2000b) and Palma et al. (2003): offset the 
RGB selection box to brighter magnitudes and count the stars 
that fall in the offset box. Assuming that the density of halo 
stars is approximately represented by a /?~ 3 power law (e.g., 
Siegel et al. 2002), the number of halo stars per unit solid 
angle is flat with respect to magnitude to first order. There- 
fore, the number of stars within the offset box should remain 
roughly constant as a function of magnitude offset. The Ca- 
rina and Ursa Minor fields were in compliance with such pre- 
dictions (Majewski et al. 2000a; Palma et al. 2003). However, 
inspection of Figure 5b already suggests that the amount of 
contamination in our case is very low because there are very 
few stars picked as giant candidates that are brighter than the 
RGB selection box. Nevertheless, we repeated the exercise 
and the results are presented in Table 2. The four faintest off- 
set bins are still counting a number of Leo I giants including 
the AGB stars left out by our color-magnitude selection box. 
At the largest magnitude offsets, the sampling of the CMD is 
incomplete due to the saturation of bright stars in our fields. 
Therefore, the bins to be considered are -5.0 < AM < -2.0. 
We take the average of these six bins and estimate the back- 
ground level in our case to be 3.5 ± 1.9 for the entire field 
and 3.4 ± 1.8 for the field outside an ellipse drawn from the 
IH95 structural parameters, which has an area of <~ 0. 1 1 deg 2 
(errors have been calculated assuming Poissonian statistics). 
This converts to a background density of 0.8 ± 0.4 deg" 2 . The 
number of Leo I giant candidates found outside the King lim- 
iting radius is twenty times this background level. While our 
derived background level is considerably smaller than was 
found in the previous studies of the Carina and Ursa Minor 
fields, Leo I is at a higher Galactic latitude and, moreover, 
we have employed a very conservative color-color selection. 
For example, if we use a different color-color selection box 
by bringing the lower-left boundary closer to the metal-rich 
dwarf locus than the one used before by A(M-D51) = 0.05, 
the background density becomes 12.5 ± 3.5 deg" 2 . Typically, 
the photometric errors are larger at fainter magnitudes. Since 
fainter Leo I RGB stars are bluer and closer to the dividing 
line between the giant and dwarf locus in the color- color di- 
agram, the likelihood of dwarfs contaminating our giant se- 
lection is greater at fainter magnitudes. Our very conservative 
selection of Leo I giant candidates gives us confidence that 
residual contamination will be very low. 

The low contamination level of our Leo I giant candidate 
list is illustrated by the upper panel of Figure 7: no Leo I giant 
candidates exist in the large area east of RA = 10! 1 22, which 
suggests that the zero-level background is reached within our 
survey region. But this begs the question of the origin of the 
numerous, very widely separated Leo I giant candidates in the 
western half of our survey area. There are 26 candidates out- 
side of twice the King limiting radius in our survey area (all 
but four of them to the west of Leo I's center), but the back- 
ground level derived above suggests that only about 3-4 of 
these stars should be contaminants. Until we can obtain spec- 
troscopy of these stars it is difficult to assess whether they are 
part of Leo I (as tidal debris), they are Milky Way contam- 
inants, or some combination of both. Tidal disruption mod- 
els (e.g. §6.1.2.) predict that stars released during the latest 
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perigalactic passage are mostly spread out along the major 
axis of the satellite, and this could account for at least half of 
these spatial outliers. Moreover, stars that became unbound 
during older perigalactic passages tend to be more irregularly 
distributed around the satellite (see Fig. 22). It is therefore 
possible that many of the Leo I giant candidates located far 
out from the main body could very well be unbound stars of 
the latter type. Certainly the large scale inhomogeneity of 
the distribution is highly suggestive of substructure of some 
kind. Nevertheless, because of the present ambiguous nature 
of these widely separated stars, we have deliberately and con- 
servatively elected to confine the remaining analyses in this 
paper to the central square degree around Leo I. 

4. STRUCTURE OF LEO I 

In this section, we present the results from profile fitting 
the two-dimensional distribution of Leo I giant candidates. A 
mass-to-light ratio of the bound population is calculated via 
the new set of structural parameters using the standard tech- 
niques of King (1966) and Ming worth (1976). We comment 
on newly identified, potential extratidal Leo I stars and also 
present isodensity contour plots of the Leo I system. 

4.1. Profile Fitting and Structural Parameters 

Structural parameters for Leo I were first derived by Hodge 
(1963, 1971), and later by IH95. All three studies used King 
profiles to fit their photographic observations. We use the 
positions on the sky of our giant candidates to explore the 
radial density profiles. To do so, we fit the surface density 
distribution of our Leo I giant candidates using two differ- 
ent models: the single-component King model and a power 
law + core (PLC). The fittings were done using a combina- 
tion of Bayesian and Maximum Likelihood techniques simi- 
lar to that employed by Kleyna et al. (1998). Errors for each 
parameter were estimated using a Bayesian approach in con- 
junction with a Markov Chain technique. The details of the 
fitting and error-estimating algorithms are described in Os- 
theimer (2002). A single-component King model has been 
widely used to fit many dSph galaxies {e.g., IH95). The PLC 
model was adopted by Kleyna et al. (1998) to fit their Ursa 
Minor surface density profile. We describe our fitting method 
below in detail. 

Several dSphs are found to have outer profile "breaks" (i.e. 
slope changes) from their central King profile density distri- 
butions: Fornax and Sculptor (Eskridge 1988a,b; Westfall et 
al. 2006), Carina (Majewski et al. 2000b, 2005; Munoz et al. 
2006), Ursa Minor (Palma et al. 2003), and Sgr (Majewski et 
al. 2003). Although statistically insignificant, IH95 found an 
excess of stars with respect to their King fit in the outer re- 
gion of Leo I. An excess of stars should show up in our radial 
density profile based on the appearance of stars spilled over 
from the main body of Leo I along the east/west in Figure 7. 
Since a fit to a distribution of stars including the excess of 
stars in the outer regions will likely inflate the derived core 
and limiting radii, we have taken the following steps to find 
a new King profile fit: First, we limit the stars used to de- 
rive the structural parameters to those within <~ 13' from the 
center of Leo I in the east/west direction. Next, six structural 
parameters (RA and Dec of the dSph center, position angle, 
core radius, ellipticity, and limiting radius) are derived from 
this stellar sample. Finally, the background level derived in 
§3.5 is combined with the six parameters to represent a newly 
fitted single-component King model. A similar approach was 
taken in the study of the Sgr profile by (Majewski et al. 2003). 



The PLC model on the other hand was fit using the entire data 
set in one pass with the background level set to the same value 
used for fitting the King model. 

Table 3 lists the best-fit parameters and errors for each pro- 
file. The first two lines in the table are for the new King and 
PLC profile fits. For comparison, we also list the King profile 
parameters derived in earlier Leo I studies by Hodge (1963), 
Hodge (1971), and IH95. The newly derived center of Leo I 
is at (a, <5)2ooo = (10 h 08 m 28 s .68, +12° 18' 19."7). 

Our derived core radius and ellipticity for either the PLC or 
King profile fit are larger than those derived by IH95, but are 
much closer to those of Hodge (1971). Similarly, our limit- 
ing radius is in better agreement with that of Hodge (1971), 
although IH95's limiting radius is also consistent within the 
errors. To illustrate how our new structural parameters better 
fit the sky distribution of Leo I giant candidates, the new King 
limiting radius is plotted over the sky distribution of Leo I gi- 
ant candidates in Figure 8 (compare to IH95 limiting radius 
size and shape in Figure 7). Figure 8 also shows the ellipse 
corresponding to the break radius (rb rea k\ see below for the 
determination of the break radius). 

With the newly derived Leo I structural parameters, we con- 
struct radial density profiles by calculating the number density 
of Leo I giant candidates in elliptical annuli of varying major 
axis radial size. Although our coverage has small gaps from 
because of spaces between the Mosaic CCD chips, our profile 
fitting program interpolates over them when doing any rele- 
vant calculation. Figure 9a shows the King model with struc- 
tural parameters derived by IH95 overlaid on the giant candi- 
date radial profile. Figures 9b show our new King model and 
Figure 9c shows the PLC profile fit. 

We now concentrate on the newly derived King model 
(Figure 9b). The observed density profile clearly departs 
from the model beyond r ma j ^10 and this "break popula- 
tion" dominates the density out to the entire survey range 
even when uncertainties in both counting and background es- 
timation are taken into account. Such profile breaks are also 
seen in N-body models of tidally disrupting satellite galax- 
ies (Johnston, Sigurdsson, & Hernquist 1999; Johnston, Choi, 
& Guhathakurta 2002). In these models, the density profiles 
of the break populations typically assume power law shapes, 
oc r" 7 . In Figure 9b, we overplot power laws with 7 = 2, 
3, 4, and 5 in the range r ma j > 10'. As seen from the figure, 
our observed profile beyond the break point is not fitted well 
by a single power law, but instead the first three radial bins 
starting at and including the break radius follow a power law 
with 7 = 3, while the next three radial bins follow a power law 
density trend with 7 = 5. Overall, a power law with 7^4 
best describes the surface density from the break radius to 
fmaj ~ 35'. It is worth noting that even for the PLC fit (Fig- 
ure 9c), the observed density profile starts to depart from the 
model at r ma j <~ 10'. This indicates that the existence of the 
break population does not strongly depend on the choice of 
fitting models. 

The surface density profile of Figure 9b roughly remains 
flat in the range 35'< r maj < 80', and then drops to the back- 
ground level beyond that. Such a flat density distribution at 
large radii has not seen in the density profiles of other dSph 
galaxies. In §6.1.1. we hypothesize that the complex form of 
the density profile might be attributable to specific patterns in 
the Leo I mass loss history. 

IH95 found that Fornax, Sculptor, and Ursa Minor 
show asymmetric residual structure after a best-fitting, 
smooth elliptical profile has been subtracted from their two- 
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dimensional distribution. Figure 8 gives some impression of 
an East-West inhomogeneity in the distribution of Leo I stars. 
To check this we plot in Figure 10 the radial density pro- 
files for the east and west halves of the galaxy, separately. In 
the inner five bins, the western half has slightly higher den- 
sities, but at low statistical significance. The densities are 
then nearly identical out to the second break. For the next 
four bins after the second break, the western half is again 
higher in density and this corresponds to the excess stars to 
the west, just outside the limiting radius. Nevertheless, the 
differences are within the Poissonian errors and we conclude 
that the structural difference between the eastern and west- 
ern halves of Leo I is of low significance. However, the fact 
that the breaks are seen independently in both east and west 
samples reinforces the reality of these structural features. 

4.2. Mass-To-Light Ratio from Core-Fitting 

If we use the King profile parameters in Table 3 to represent 
the bound, part of Leo I, the bound mass can be estimated un- 
der the usual assumptions of virial equilibrium and the stan- 
dard equation given by Ming worth (1976): 



where R Ctg is the geometric mean of the King core radius 

(= r c y/l — e) in parsecs (271 ± 26 pc for Leo I), fx is the King 
(1966) mass parameter, and j3 is a velocity parameter which 
is related to the observed velocity dispersion. Both fi and 
(3 are strongly dependent on the concentration of a system 
(King 1962). From the results in Table 3, the concentration of 
Leo I is log(r, /r c ) = 0.39 ± 0.04. The scaling parameter /i was 
taken from an extrapolation of Table II of King (1966) and [3 
from Figure 4-11 of Binney & Tremaine (1987). The cen- 
tral velocity dispersion of ao = 8.8 ± 1.3 km/s (M98) yields 
M tot = 4.0(±1.2) x 10 7 M Q 16 . To obtain the corresponding lu- 
minosity, we integrate our best-fit model King profile for giant 
star candidates and scale it to the total luminosity using a cen- 
tral surface brightness of /io,v = 22.4 mag/arcsec 2 taken from 
Table4 of M98; this yields L tot = 7.6 x 1O 6 L . The M/L using 
this technique thus becomes (M/L), ot y = 5.3 ± 1.6M Q /L Qi y, 
where the error is from the uncertainty in the mass derived 
above. Though we have used different structural parameters, 
we obtain nearly the same M/L as M98. We note that this 
M/L is actually not that large compared to those typically 
found for dSphs, and is fairly similar to that of dE galaxies 
and globular clusters of similar luminosity. As has been pre- 
viously found (M98), Leo I has one of the lowest M/L among 
the Galactic satellites, though one that, by this method, still 
suggests a significant amount of DM. 

4.3. Isodensity Contours 

In Figure 11, we present the isodensity contour map of 
Leo I. This map was constructed using the Leo I giant can- 
didates in the following manner: The equatorial coordinates 
of each star were first converted to a flat Cartesian system 
via tangential projection centered at the newly derived system 
center. The Cartesian space was divided into a large grid at 
intervals of 2.'7 in each dimension. We then counted the num- 
ber of stars in each grid interval to obtain the density for that 
part of the map. The gaps between the chips of Mosaic CCD 

16 Our own derived central velocity dispersion of 8.2 km s -1 (see §5.5) 
only lowers the estimated Leo I mass by 14%. 



camera are 50 CCD pixels in the rows and 35 CCD pixels in 
the columns, corresponding to <~ 13" and ~ 9" in the sky, re- 
spectively. In order to check whether the absence of stars in 
the CCD gaps produce any artifacts in Figure 11a, we per- 
formed checks on our data as follows: First, we divided the 
CCD gaps into segments of 6' bins, and filled each bin with 
randomly placed artificial stars. The number of stars that went 
into each bin was calculated by taking the average of the num- 
ber of Leo I stars on either side of the gaps in similar RA/Dec 
ranges. Once the gaps were filled, an isodensity contour plot 
was constructed using the exact same procedure as described 
above. We have also placed artificial gaps with same size as 
the real ones in random places and constructed several contour 
plots. Only slight differences were found between the isoden- 
sity contours constructed in these ways and Figure 11a. This 
ensures us that the influence of the CCD gaps on Figure 11a 
is negligible. 

In order to reduce the granularity from the griding process 
in our isodensity contour, we constructed a grid four times 
finer than the one used above, and replaced each star with a 
6.0 x 6.0 arcmin 2 grid filled with a two-dimensional Gaus- 
sian template (a = 1.1 arcmin; amplitude = 1). We assigned 
density to each pixel by counting up the numbers in a simi- 
lar manner discussed above. The resulting smoothed image is 
shown in Figure lib. 

Overall, the inner contours of Leo I are similar to the iso- 
pleth map shown in IH95 (see their Fig. Id). However, the 
spill-over of stars at the profile break radius is evident in the 
outer contours of Figure 1 1 as contour extensions stretching 
out from the main body along the semi-major axis. It is impor- 
tant to note that the spatial density for the outermost contour 
in Figure 11a is still higher than 600 times the background 
level derived in Section 3.5. 

5. SPECTROSCOPY OF LEO I GIANT CANDIDATES 
5.1. Observations and Basic Reductions 

We obtained spectra of a total of 135 stars in the Leo I field; 
49 stars lying just inside and outside the King limiting radius 
to the west of the Leo I center (Field KD1), and 86 stars near 
the center and along the east major axis of Leo I (Field KD2). 
The goals of our spectroscopic observations are (1) to check 
further the reliability of our color-color and color-magnitude 
selections discussed in §3, (2) to verify the reality of the break 
population seen in the radial density profiles of Leo I (see 
Figure 9), and (3) to study the dynamics of both the inner and 
outer parts of the Leo I dwarf. 

The spectra were obtained with observations of two multi- 
slit masks with the Deep Imaging Multi-Object Spectrograph 
(DEIMOS; Faber et al. 2003) on the nights of UT 2003 Oc- 
tober 29 (Field KD1) and 2004 October 15-18 (Field KD2) 
using the Keck II 10 m telescope (see Figure 16 for the place- 
ment of our masks on the sky with respect to the Leo I giant 
candidates). The masks were designed using coordinates from 
our photometric catalog, which have been locked to the astro- 
metric system of the USNO-A V2.0 catalog (see §2). Priority 
for slit selection was given to stars selected to be Leo I giant 
candidates by the criteria described in §3. Even for this prior- 
ity sample, all stars falling within the 5'x 16' slitmask area of 
Field KD1 could not be observed because of slit overlap. On 
the other hand, occasional gaps in mask area remained and 
these were filled with slits for, in priority order: (1) stars se- 
lected as giants, but not as Leo I giants, and (2) any other star 
in the field without regard to classification. Observations of 
stars classified as dwarfs prove useful as checks on the reli- 
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ability of the dwarf/giant separation. The spectrograph was 
configured with the 1200 lines/mm grating and 1 "slits with 
the central wavelength set at 7800A. This instrumental set up 
provides 0.33A/pixel dispersion, a spectral resolution FWHM 
of 1 .95 A (68 km s -1 at the Ca infrared triplet) after accounting 
for the 0.7 x anamorphic magnification factor, and a spectral 
coverage of 6500-9 100 A (which varies a little from slit to 
slit depending on their particular placement within the mask 
field of view). The total integration times were 80 minutes 
for Field KD1 and 115 minutes for Field KD2, divided into 
four and five separate exposures, respectively, for cosmic ray 
removal. Typical DEIMOS spectra for different magnitude 
stars are shown in Figure 12. The typical seeing during the 
Leo I spectroscopic observations were 0."8 - 1 ."0. All steps of 
the data reduction to wavelength-calibrated, one-dimensional 
spectra (including bias subtraction, flat fielding, sky subtrac- 
tion, wavelength solutions, and extraction of 1 -d spectra) were 
done using the spec2d reduction pipeline (Cooper et al., in 
prep) developed at the University of California at Berkeley 
for the DEEP2 galaxy redshift survey project. 

5.2. Radial Velocity Determination 

Derivation of the radial velocity (RV) measurements uses 
a sequence of steps, each devised to address some aspect of 
noise suppression or bias compensation. Combined, the re- 
sulting scatter in the achieved RV measurements can be shown 
to be less than 1/40 of the FWHM of a resolution element for 
the best exposed Leo I stars. 

The first phase of the RV reduction follows the masked, 
Fourier-filtered cross-correlation process outlined in Majew- 
ski et al. (2004) and used in a number of our RV studies over 
more than a decade (e.g., Kunkel, Irwin, & Demers 1997). 
Briefly, each star is first Fourier-filtered to attenuate frequency 
components lower than those given by the typical absorption 
line and higher than permitted by the intrinsic line resolution 
(as might be introduced by cosmic rays, for example). Mask- 
ing multiplies a filtered spectrum in wavelength space with 
zero to suppress, or with unity to admit portions of a spectrum. 
Mask wavelengths are selected according to several line lists 
that each serve specific functions. For example, an OH night 
sky emission mask multiplies all candidate data spectra prior 
to the correlation process to suppress possible residual noise 
due to imperfect night sky subtraction. Another mask simi- 
larly multiplies filtered spectra prior to correlation to suppress 
telluric water and O2 bands when their strength at the instru- 
mental resolution exceeds a few percent. Finally, a third line 
list is used to mask the cross-correlation "master template", 
which is set to zero everywhere except at low-ionization or 
low-excitation metallic or Balmer features that, at the instru- 
mental resolution and at mean target star signal levels, con- 
tribute a useful equivalent width component for a metallic - 
ity representative of the anticipated population (in this case, 
Leo I). The suitability of features selected for inclusion in this 
last mask line list is decided from visual inspection of the Arc- 
turus Atlas (Hinkle et al. 2000), which permits identification 
of satisfactory features as well as untrustworthy blends (usu- 
ally from line pairs of dissimilar elements) for exclusion, and 
with reference to representative target star spectra to deter- 
mine in a practical sense the visibility of the selected features. 

In the specific case of our Leo I reductions, our master 
cross-correlation template was based on the star 61788 in the 
KD1 mask, which was very well-exposed (~ 10,000 ADU 
pixel" 1 — roughly twenty times stronger than the weaker 



Leo I spectra; see Fig 12). This source was used as an RV tem- 
plate because no exposure of a more standard RV reference 
star was obtained, since the Keck observations of Leo I were 
obtained only as a backup program for the DEEP project. The 
benefit of using this particular star as a cross-correlation tem- 
plate is its close match in spectral type and metallicity with 
the Leo I spectra. The disadvantage is that the true RV of this 
local reference is unknown, so the final RVs derived in the 
initial cross-correlation phase are subject to zero-point uncer- 
tainties, although their random errors are low. Since the key 
dynamical descriptors of the Leo I population of relevance to 
our later analysis are based on relative distributions of RVs, 
any systemic RV offset present in the data does not explic- 
itly affect our interpretations. The correction of systematic 
velocity offsets to a known velocity system can be achieved 
by accounting for offsets against telluric line references, as 
described below. In addition, we describe how we ultimately 
adjust to the zero-point of the M98 RV system in §5.3, a nec- 
essary additional step if we wish to combine the two data sets. 

In preparing star 61788 to serve as a cross-correlation tem- 
plate, the masking procedure left predominantly the lines of 
the Ca infrared triplet, as well as the stronger features of Mg, 
Fe, Ti, Ni, Si, Na, and Ha. These lines visible in the atmo- 
sphere of the star 61788 were preserved with the non-zero 
mask "slats" represented by a Gaussian of unit amplitude and 
a width roughly 1.2 times the instrument profile. A total of 
85 lines were used in cross-correlations using this masked 
template, with each slat having a FWHM of 3. OA for a to- 
tal used spectral range of 255A over the full spectral range. 
Experiments using slats of 2.5A and 4.2A yielded no signif- 
icant variations in the quality of the cross-correlation results. 
As described above, before cross-correlation the other Leo I 
stars have been similarly Fourier-filtered to the 61788 spec- 
trum. To describe the quality of each cross-correlation, we 
have adopted the RV quality index Q, which is a descriptor 
of the shape of the central cross-correlation peak with respect 
to sideband peaks. 17 We note that although the meaning of 
the RV quality values are identical to those given in Majewski 
et al. (2004), the actual values of the cross-correlation peak 
(CCP) levels are on a relative scale that is only meaningful 
within the context of the present instrument set up. In gen- 
eral, only those RVs with Q>4 can be trusted as reliable and 
we discard any star that does not meet this criterion in subse- 
quent analysis. 

A primary limitation on the precision of RVs as just ob- 
tained is the non-uniform manner in which a spectrometer slit 
is illuminated, even when the star is perfectly centered in the 
slit. The presumption that a spectrometer slit is uniformly il- 
luminated is true only the when illuminating source is uni- 
formly spread over the portion of sky that the telescope is 
sampling. Even in poor seeing a stellar point source, through 
spread over a angle with a roughly Gaussian profile on the sky 
tends to illuminate the slit aperture (along dispersion) differ- 
ently compared to the illumination from a lamp used to im- 
pose pieces of monochromatic light for wavelength calibra- 
tion. To the extent that this difference appears perpendicular 
to the slit aperture, there is a displacement in the centroid of 
(monochromatic) stellar light along the dispersion compared 
with that of light from a uniformly illuminated source (such 
as the sky or a calibrating lamp). The effects are obviously 
greatly amplified in the case that the star is not centered in 

17 See Majewski et al. (2004) and Kunkel, Irwin, & Demers (1997) for 
descriptions on how Q values are assigned. 
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the slit, whether due to random errors in the astrometry, im- 
perfect tooling of the mask slits, and/or mask misalignment at 
the time of observation. The dispersion displacement varies 
with the circumstances of an exposure and can be determined 
only from signatures carried by the stellar signal itself. Fortu- 
nately, the signature of this displacement is also imposed on 
absorption features originating from molecules of telluric O2 
and H2O in the light path, and whose profiles and positions ac- 
curately track the mean integrated slit function of each partic- 
ular star as passed through its slit. Thus, the mean wavelength 
offsets in telluric features can be used to derive slit displace- 
ments on a star-by-star basis. To derive the telluric offsets 
we design a new cross-correlation mask in which all spectral 
features are blocked except the regions with the strongest tel- 
luric absorption features. The wavelength ranges upon which 
the telluric offset determination is based are 6866-69 12A (the 
Fraunhofer B band), 7167-7320A, 7593-7690A (the Fraun- 
hofer A band), 81 10-8320A, and 8925-9120A, but all telluric 
features weaker than 20% were not used. A "master" tem- 
plate of telluric features was made from a fairly strong Keck 
spectrum of a star with weak stellar features. Unlike stellar 
features, the telluric features will look identical in all stars 
(modulo the S/N of their spectra), and so the dominant peak 
profile in cross-correlations with our masked telluric spectrum 
should appear the same. This was confirmed visually, and 
the lowest "Q" values we would assign to any of these cross- 
correlation peaks would be a "6"; even stellar spectra with 
bad Q from the cross-correlation of stellar features still pro- 
vide excellent telluric correlation profiles. We note that appli- 
cation of the offsets derived from cross-correlation of telluric 
features also corrects the original RVs for errors in the pixel- 
to-wavelength calibration and, in addition, places the RVs to 
an absolute reference. 

The heliocentric RVs corrected for these offsets (with the 
additional zero-point correction to the Leo I systemic velocity 
as derived in §5.3) are given in column (8) of Table 4, along 
with their RA and Dec coordinates, Mo magnitudes, (M— Ti)o 
and (M-D51)o colors in columns (3) through (7). We also 
list Galactocentric standard of rest (GSR) RVs 18 , telluric off- 
sets, cross-correlation peaks (CCP), and quality index (Q) in 
columns (9) through (12) of Table 4. 

For those stars observed in 1" wide slits 19 , the mean tel- 
luric offset is -1 1 .5 km s" 1 with a dispersion of 7.7 km s" 1 for 
those stars in mask KD1 and -2.8 km s" 1 with a dispersion of 
3.5 km s" 1 for those stars in mask KD2. That both mean off- 
sets are non-zero is a reflection of systematics in the original 
template derived from star 61788; that the mean offsets are 
different for each mask reflects the differences in position- 
ing of each mask relative to their respective star fields. The 
7.7 and 3.5 km s" 1 dispersions are respectively equivalent to 
~ 0.2 and ~ 0.1 arcsec variations in the slit centering, and re- 
flects the quality of the original astrometry (see §2), the slit 
manufacture and the rotational alignment of the slit mask on 
the sky. Figure 13, which shows the telluric offset correction 
for both masks as a function of position along the mask (ap- 
proximately declination for mask KD 1 and right ascension for 

18 We convert to Galactocentric standard of rest by assuming a local 
standard of rest velocity of 220 km s -1 and a solar peculiar velocity of 
(M,v,w) = (-9,12,7)kms- 1 . 

19 Seven stars listed in Table 4 were actually alignment stars observed in 
4"wide slits. However, with the telluric offset corrections we can actually 
derive good RVs for these spectra with larger equivalent slit functions. The 
larger telluric offset corrections needed for these stars are evident in the tab- 
ulated values. 



mask KD2), demonstrates however that the initial astromet- 
ric reduction contributes significantly to the star-by-star slit 
miscentering, given the clear correlations of the offset trends 
with the original CCD frame on which the star was originally 
photometered (each CCD frame has a unique astrometric so- 
lution). The dispersions in the telluric offsets represent the 
actual precisions limits for multi-slit RVs in the case when 
the telluric offsets are not accounted for. After applying the 
telluric offset the velocity precision is therefore well better 
than this. 

5.3. Determination of Sampling Errors 

Based on the relative strengths of the cross-correlations 
against the stellar and telluric templates, the dominant con- 
tribution to the RV errors are in the stellar absorption cross- 
correlation; the telluric cross-correlation peaks are almost al- 
ways much stronger than the kinematic correlation peaks, 
because the equivalent width of telluric absorption features 
sum to an order of magnitude larger than those intrinsic to 
the stellar atmosphere. Errors in the stellar absorption cross- 
correlation are partly a reflection of differences in stellar line 
strengths between the template and target spectra. 

An assessment of random sampling errors may be ap- 
proached by a variety of techniques; however, the method of 
Tonry & Davis (1979) is not one of them. The adopted cross- 
correlation methodology employed here differs from that of 
Tonry & Davis, having evolved over years of experience and 
motivated by the goal to cover as broad a range of (generally 
unknown) spectral types among the target stars as possible. 
The Tonry & Davis method relies explicitly on the character- 
istic that a template and a target spectrum must be of identical 
shape, or spectral type, so that the resulting cross-correlation 
ideally becomes an even function with null imaginary terms. 
Then the imaginary terms that appear for real data may be uti- 
lized for the estimation of sampling errors. In contrast, our 
methodology, by using slat masking, maximizes the tolerance 
of the cross-correlation to widely variable spectral types, so 
that F through M stars have in the past been measured equally 
successfully to full precision in the same correlation process 
with a single template design. But such dissimilarities invali- 
date the reliance on the even character of the correlation func- 
tion, and imaginary terms are generally non-zero even when 
sampling errors are rigorously zero. Consequently alternate 
methods for estimating errors with our methodology are re- 
quired. 

To obtain a first estimate of the errors in our RVs, we 
repeated the stellar absorption cross-correlations twice with 
other RV templates manufactured from other stars in the Leo I 
masks with much stronger stellar absorption features than 
stars 61788, namely stars 61736, 61782. With the cross- 
correlations there was little difference in the final results for 
the best quality target spectra, and there were only minor 
changes in Q values using different RV templates. For Q = l 
spectra the RMS scatter in RVs was 1.2 km s -1 . This rises 
to 3.0 km s" 1 when any one of the cross-correlation peaks 
dropped to Q = 4. This analysis suggested an uncertainty scale 
of 1.2, 1.5, 2.0, and 3.0 km s" 1 for stars with Q = 1, 6, 5, 
and 4, respectively. However, these uncertainties must be re- 
garded as lower limits, since the use of different templates on 
the same target spectra do not lead to completely statistically 
independent measurements. 

Because repeat observations of our target stars are not avail- 
able, an alternative error estimation derives from creating ar- 
tificial "repeat observations" via the convolution of the raw, 
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unfiltered candidate spectra with randomized Poisson noise. 
Practically, this is accomplished by using the original spec- 
trum as a probability distribution function, whereby, at any 
wavelength, the envelope of the original spectrum defines 
the "expectation" of received source counts with respect to 
neighboring wavelengths. With each wavelength "bin" set 
to the resolution FWHM width, the spectrum can be repopu- 
lated with Poisson distributed photon events, dropped in ran- 
domly with wavelength destinations weighted according to 
the expectation probabilities. To get a representative spectrum 
whose noise characteristics are like those of the raw data (pre- 
suming a negligible contribution from read noise) one keeps 
adding Poisson events until the count of events equals that of 
the raw target spectrum. On average, each wavelength bin in 
the new pseudo-spectrum retains the same amplitude, com- 
pared to neighbors, as the original spectrum, with modulation 
only due to Poisson noise. When a number of pseudo-spectra 
are so rebuilt but arbitrarily scaled to a much stronger total 
photon count than the original spectrum signal strength, their 
cross-correlations with the template show zero velocity scat- 
ter, as one expects for strong signals. When, on the other 
hand, Poisson distributed photons are dropped into each slat 
so that the overall mean exposure level of the pseudo-spectra 
are comparable to the observed spectrum, the scatter in the 
derived velocities is found to be identical to that obtained by 
re-observing the same star as many times. 20 

The above error estimation method was applied to a vari- 
ety of Leo I spectra with the generated pseudo-spectra cross- 
correlated against the same RV template using the same 
masks. This procedure packages into one procedure a test 
of all aspects of the data processing stream, including all in- 
strumental traits, as well as the design of the line-lists, tem- 
plate fabrication, etc. If we adopt test spectra near the low end 
limit of signal levels, i.e. <~ 500 ADU, we obtain dispersions 
of 1.55, 1.77, 2.87 and 4.91 km s" 1 for stellar feature cross- 
correlations in Q= 7, 6, 5 and 4 spectra, respectively. Obvi- 
ously, better precisions are achieved for the large number of 
spectra with higher signal levels; e.g., for Q = l spectra, dis- 
persions of 0.44, 0.50, 1.46 km s" 1 were obtained for signal 
levels of about 4200, 2500, and 1000 ADU. Independent tests 
of the telluric cross-correlations at these signal strengths lead 
to scatters of about 1.2 km s" 1 (and significantly better for 
stronger spectra). Adding these contributions in quadrature, 
and conservatively adopting the worst case (i.e. low signal 
strength) dispersions for each quality, leads to the following 
estimated uncertainties for stars with Q = 7, 6, 5 and 4, respec- 
tively: 2.0, 2.1, 3.1, and 5.1 km s -1 . We adopt these as our 
velocity uncertainties for subsequent analyses. However, it is 
worth pointing out that by adopting conservative estimates of 
our uncertainties at every level, we do risk underestimating 
derived true velocity dispersions after we remove the contri- 
butions from observational uncertainties. 21 

5.4. Testing the Photometric Selection of Leo I Giants 

20 The basis for this claim derives from tests of this methodology on other 
data sets available to the authors wherein multiple observations of the stars are 
available. In all cases tested, equivalent RV scatters were found for multiple 
observations of the same star compared to multiple pseudo-spectra generated 
from one observation. 

21 For example, adopting the lower limits to the velocity uncertainties de- 
rived from our first estimate above (i.e. 1.2, 1.5, 2.0 and 3.0 km s _1 for Q = 7 
to 4) raises the central Leo I velocity dispersion (see §5) by about 5% (from 
8.2 to 8.6 km s~') and the "cold point" velocity dispersion near the King 
limiting radius (see §5) by the same fraction (from 3.9 to 4.1 km s~'). 



The extreme systemic RV of Leo I, v/, e ; = 287 km s" 1 (M98) 
is a great advantage for distinguishing true Leo I members 
from Galactic stars — very few Galactic stars are expected 
to have such extreme velocities at this position in the Galaxy 
[(l,b) = (226, +49)°]. Figure 14 shows a histogram of the ob- 
served RVs of stars listed in Table 4; stars photometrically se- 
lected to be Leo I giant candidates are shaded darkly. Leo I gi- 
ants are easily identified by the clump at Vhd / (km s" 1 ) > 250. 
In KD1, 21 out of 42 stars are identified as Leo I members and 
in KD2 all 83 stars are Leo I members, based on their RVs. 
Stars we targeted simply as "mask fillers" tend to have RVs 
clustering near km s" 1 as expected for MW disk stars. 

In Figure 15, we plot the stars observed with the Keck 
DEIMOS over the color-color and color-magnitude diagrams 
with the selection regions shown. In both the color-magnitude 
and color-color diagrams, the stars with measured velocity 
Vhei > 250 km s" 1 are generally located within our primary 
giant selection region whereas the stars with lower velocities 
are not; in particular, in the two-color diagram the lower ve- 
locity stars all fall along the elbow-shaped, dwarf star locus, 
suggesting these stars to all be relatively nearby Milky Way 
dwarfs. Even more gratifying is that for those 96 stars se- 
lected to be Leo I giants using both the color-color and color- 
magnitude diagram and for which we obtained RVs, all of 
them are Leo I members by the RV criterion; this suggests 
that the reliability of our selection is very good, and we em- 
phasize that a significant portion of our sample includes stars 
selected in the very low density outer parts of the Leo I system 
(see Figure 16). We also observe that our Leo I giant can- 
didate selection criteria are on the conservative side, because 
nine additional stars that we did not pick as giant candidates 
are also Leo I RV members. In most cases a slight expansion 
of our selection criteria would have picked up all but one 22 of 
these other Leo I stars, but this expansion would have brought 
in false positives. In particular, it may be seen that our catalog 
can be trusted to even fainter magnitudes than the conserva- 
tive Mo < 21.5 limit we adopted for the analyses presented 
in our photometric studies. It may also be seen that use of 
the M-D5 1 color eliminates contaminant dwarf stars that fall 
within the color-magnitude locus of the Leo I RGB, and thus 
improves the efficiency of our target selection. This reliabil- 
ity analysis lends credibility to the structural features we have 
mapped using Leo I giant candidates selected in this way. 

Figure 16 shows the sky distribution of the stars in Fig- 
ure 15. Out of the total 105 Leo I giants observed with Keck 
DEIMOS, 90 stars lie within the King limiting radius and 48 
stars lie within the core radius, leaving 15 confirmed Leo I 
members beyond the limiting radius - i.e., in the domain of 
the "break" population we identified in our study of the Leo I 
light profile. The new Keck spectroscopy confirms that our se- 
lection of Leo I giant candidates is both effective and reliable, 
and so lends confidence that the population of stars falling 
outside the King limiting radius (which may be the nominal 
tidal radius of the system - see §6.1 below) and constitut- 
ing our detected "break population" (e.g., Fig. 9) are actually 
stars associated with Leo I. 

5.5. Velocity Dispersion 

The variation of velocity dispersion (<j v ) with projected ra- 
dius of dSphs is important for testing dynamical models of 
dSphs (e.g., Kroupa 1997; Lokas, Mamon, & Prada 2005; 

22 This one Leo I RV member well outside our color-magnitude selection 
box at (M-T 2 )o < 1.0 (see Figure 15) is likely to be an anomalous Cepheid. 
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Strigari et al. 2007; Munoz, Majewski, & Johnston 2007). 
For a bound stellar system where mass-follows-light <j v is ex- 
pected to decrease with radius and approach zero at the tidal 
radius. On the other hand, dSph models with extended DM 
halos predict a dispersion that falls off more slowly than mass- 
follows-light models (e.g., Kleyna et al. 2001; Wilkinson et 
al. 2002). Finally, tidal disruption models predict fiat/rising 
<r v profiles (e.g., Kroupa 1997; Munoz et al., in preparation) 
into the domain dominated by stars that have become unbound 
during/after tidal interactions with the host galaxy. 

With our Keck DEIMOS data, we can explore the radial 
variation of <t v out to and beyond the King limiting radius. 
Before doing so, we adjust our RVs to the zero point of the 
M98 RV system so that both data sets may be considered to- 
gether for improved statistics. The weighted mean a v of the 
33 stars observed by M98 is 8.6 ± 1.2 km s" 1 . These stars 
lie in the semi-major axis range r ma j/ru m = 0.055-0.367 ac- 
cording to our newly derived structural parameters. The a v 
for 50 stars we observed in the same r maJ range is 8.2 ±0.9 
km s" 1 , which is consistent with the M98's velocity dispersion 
within the error. 23 We take as the offset to our RV system the 
+12.4 km s" 1 difference between the error-weighted mean of 
our and the M98 samples, and apply this offset to the entire 
Keck DEIMOS data set (although it is not clear, in fact, which 
of the two surveys is closer to the true absolute RV system). 
In the upper panel of Figure 17, we plot the heliocentric RVs 
versus the projected major-axis radial distances for both Keck 
DEIMOS and M98 samples. For comparison, we have plotted 
each sample with different symbols (and colors). The remark- 
able similarity of the RV distributions within 5 arcmin from 
the center of Leo I assures us that both samples are on the 
same system. We also show the mean RVs along projected 
radius in the lower panel of Figure 17. The mean RV trend 
shows that there is a hint of rotation in the inner < 5 arcmin, 
but with low statistical significance. Spectroscopic observa- 
tions for stars to the west in the range r ma j = 5-10 arcmin 
will help reveal whether rotation is significant for Leo I. We 
use the combined sample of M98 and Keck spectroscopy in 
subsequent analyses. 

In Figure 18, we present Vhd and a v as a function of ra- 
dial distance from the Leo I center, calculated for elliptical 
(left panels) and circular radii (right panels). We use differ- 
ent number of stars per bin for the middle and lower panels. 
The overall trend of the a v is an initial decline followed by a 
flat or rising profile. Similar profile behaviors up to and be- 
yond the nominal King radius are also seen in the velocity 
dispersion profiles of Ursa Minor (Munoz et al. 2005), Draco 
(Munoz et al. 2005), and Sculptor (Westfall et al. 2006). As 
discussed in these references, such behavior means either that 
we are seeing evidence for tidal disruption in these dSphs, or 
that these satellites have very extended DM components dy- 
namically traced by these stars at large radii. 

Closer inspection of the upper panels of Figure 18 shows a 
number of stars (more at larger radius) that have an RV de- 
viating significantly from the mean. These outliers 24 are pre- 
dominantly toward the high RV side, and create an obvious 
asymmetric tail in the distribution of Leo I velocities shown 
in Figure 14. Indeed, as shown in Figure 19, it appears that 
these outliers create a secondary hump in the RV distribution 

23 All velocity dispersion calculations throughout this paper follow the 
method of Armandroff & Da Costa (1986). 

24 Although we use the term outliers, these are obviously associated with 
Leo I. 



about 20 km s" 1 higher than the Leo I mean RV of 288.8 km 
s" 1 , but in the least there is an obvious asymmetric tail towards 
larger RVs. 

Members of the outlier population we identify as those stars 
with RVs beyond 2.3-ct (where sigma refers to the velocity 
dispersions of the entire sample, i.e. 8.9 km s" 1 ) for the Leo I 
RV sample with r ma j/ru m < 0.6, and those with RVs beyond 
1.5-cr for r ma j/ru m > 0.6 Leo I stars. Our tighter criterion 
for the stars at larger radii is based on examination of the 
RV distribution shown in Figure 19: the tight peak in the 
rmaj/rum > 0.6 distribution clearly stands out from a broader 
distribution of RVs and is well approximated by a Gaussian of 
dispersion of 3.9 km s~ l . Because this tight peak is centered 
on the mean velocity of the Leo I core, it seems appropri- 
ate to associate it with the typical RV population at smaller 
radii (though the latter has a broader distribution). But it is 
apparent that about half of the outer Leo I stars are part of a 
second, much more broadly RV-distributed, "outlier" popula- 
tion in Figure 19b. Adopting a 3-er cut based on the rather 
cold, 3.9 km s -1 central dispersion peak among the stars with 
r m aj/riim > 0.6 identifies the same set of RV outliers as a 1.5-a 
scheme where the entire Leo I sample a is used. 

The outliers so identified are marked with open circles in 
the top panels of Figure 18. The spatial distribution of these 
stars is shown in Figure 20 with respect to the entire sample 
of Leo I giant candidates, where it may be seen that not only 
do they tend to lie to the west side of Leo I, they are predom- 
inantly found near or beyond the nominal King radius. The 
latter point is not simply due to the combined facts that we 
have used a tighter "outlier" definition for stars at large radii 
and that all of the stars in the western DEIMOS mask are at 
large radii: In fact, there are still 1 1 stars to the east side of 
Leo I center with r ma j/ru m > 0.6, but only one of these is an 
outlier by the above definition, compared to 9 outliers among 
22 r ma j/ru m > 0.6 stars on the west side. Of the 9 outliers to 
the west, eight are outliers to the positive RV side. 

Exclusion of these outliers from the calculation of a v trends 
leads to different results, as shown by open circles in the mid- 
dle and lower panels of Figure 18: the velocity dispersion 
tends to fall with radius, and to a relatively cold value at large 
radius. For the stars beyond the King limiting radius, the dy- 
namically cold population has a mean RV of 287.7 ± 1 .6 km 
s -1 , close to the systemic velocity of the Leo I core, and a ve- 
locity dispersion of 4.0 ± 1 .6 km s -1 . These dynamically cold 
stars also appear to be more in line with the Leo I major axis, 
as seen in Figure 20. The RV outlier stars, when included, in- 
crease the dispersion outside the limiting radius to 10.3 ± 2.2 
km s . 

Both the presence and meaning of velocity dispersion "cold 
points" near the nominal King radii of various dSphs has been 
debated recently (Wilkinson et al. 2004; Lokas, Mamon, & 
Prada 2005; Munoz et al. 2005). The reality of the phe- 
nomenon turns on how one deals with (i.e., includes or ex- 
cludes) apparent outliers and how one bins the data (Lokas, 
Mamon, & Prada 2005; Munoz et al. 2005). In the present 
case, where the Leo I systemic velocity is so extreme, it 
is hard to believe that the outliers are anything but Leo I- 
associated stars, and so we are presented with a situation 
where the outliers and the dispersion characteristics demand 
an explanation within the context of a dSph structural model. 

The lopsided nature of the RV outliers (i.e., that they tend 
to be outliers to the high RV side) provides a significant clue 
to the nature of the outer parts of the Leo I dSph. Were the RV 
outliers bound stars within a large Leo I DM halo, it is hard to 
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understand why they would not exhibit symmetrical dynamics 
— i.e., a hot population of stars with members both approach- 
ing and receding relative to the Leo I core. The lack of any 
particular spatial concentration of the RV outliers (other than 
tending to the west side of Leo I) does not support the idea that 
the RV asymmetry is caused by a star cluster or "dynamical 
fossil" like that recently suggested to be "sloshing back and 
forth within the" DM halo of the Ursa Minor dSph (Kleyna 
et al. 2003). The observation of significant spatial "substruc- 
ture" within dSphs (e.g., Olszewski & Aaronson 1985; De- 
mers et al. 1995; Eskridge & Schweitzer 2001; Kleyna et al. 
2003; Palma et al. 2003; Coleman et al. 2004; Coleman & 
Da Costa 2005; Olszewski et al. 2005) has long been a con- 
cern for DM models, where high DM densities should quickly 
smooth such substructure out. Here we are presented with an 
apparent substructure evident both dynamically and spatially. 

A more natural explanation of the RV distribution is that the 
RV outliers and many of the stars outside the nominal King 
limiting radius represent stars that have likely been tidally 
stripped from Leo I, whereas the more dynamically cold com- 
ponent at the nominal Leo I systemic velocity seen at large 
radii represents the outermost reaches of the bound population 
of Leo I stars. The decreasing RV dispersion trend when the 
outliers are excluded suggests that we are seeing the tapering 
gravitational field of a Leo I apparently lacking an extended 
DM halo. Moreover, the observed asymmetric RV distribu- 
tion is naturally produced by tidal debris, as we show in §6.3. 
A spectroscopic survey over a larger area would of course help 
clarify the dynamics of the outer Leo I system — in particu- 
lar, spectroscopic observations at large radius on the eastern 
side of Leo I would provide particularly strong leverage on the 
tidal debris model, which would predict negative RV outliers 
to the east side of Leo I. 

6. DISCUSSION 

6.1. Dark Matter, Extended Halo s and Tidal Stripping 

The question of the origin of the apparently large M/L in 
dSph galaxies remains one of the most vexing in Local Group 
research. While the presence of large amounts of DM is the 
most popular explanation (e.g., M98), other proposals, includ- 
ing Modified Newtonian Dynamics (MOND; see Sanders & 
McGaugh 2002, and references therein), tidal heating and/or 
disruption, dynamical resonances (Kuhn & Miller 1989), and 
even the notion that dSphs are completely unbound (Kroupa 
1997; Klessen & Kroupa 1998) have been proposed to ex- 
plain the large velocity dispersions observed in the cores of 
the Galactic dSphs. Leo I provides an interesting opportunity 
to revisit some of these mechanisms, because its large Galac- 
tocentric distance and velocity implies that Leo I could only 
have had a few encounters with the MW center. Tidal heat- 
ing/disruption or resonance effects, if they occur, will have 
been limited to those that could organize and be sustained 
over only one or two perigalactica. Because of this, Leo I 
offers the unique chance to gauge the importance of such ef- 
fects on a per perigalacticon basis. 

6.1.1. The Case for Tidal Disruption 

In §4.2 we have rederived the M/L of Leo I using stan- 
dard core-fitting prescriptions and the central velocity disper- 
sion. The results are not much different than previous find- 
ings, which suggest a modest DM component in Leo I com- 
pared to that in other dSphs. A primary point of this paper 
is to establish whether, independent of the net DM content of 



the dSph, it may have extended structure indicating that it is 
being tidally disrupted. 

We believe the weight of evidence suggests compellingly 
that Leo I has indeed been tidally disrupted and, moreover, 
that the derived King limiting radius approximates the true 
tidal radius in this particular system: 

( 1 ) We find a significant number of widely placed giant can- 
didates associated with Leo I in our survey area (see Figure 
7), and with an especially pronounced density near the King 
limiting radius. While it might be argued that these stars rep- 
resent a secondary, bound "halo" population around Leo I — 
perhaps tracing an extended DM halo — as can be seen in Fig- 
ure 8 the "break population" stars are more spread out along 
the east-west direction than in north-south in a manner that 
resembles tidal arms. In order to demonstrate the "bipolarity" 
of these outer stars more clearly, we divide the region out- 
side the break radius of Leo I into azimuthal sectors of 18° in 
width and 17' in major axis radial length (r ma j = 10f4-27.'5), 
and count the giant candidates in each sector 25 . The counts 
were then normalized by the surveyed area in each sector to 
obtain the mean densities as a function of azimuthal sector 
presented in Figure 21. As may be seen the "arms" are repre- 
sented by two broad density peaks separated by ~ 180° to the 
east and west. A minimum in density is apparent to the north, 
but a corresponding minimum to the south is interrupted by 
the presence of a small peak at 145-180° which originates 
from the bridge of seven stars extending to the southeast (vis- 
ible in Fig. 8). Whether this feature is a chance alignment of 
Leo I stars or a real dynamical structure may require high res- 
olution velocity data to resolve, but we note that if these stars 
were a bit more spread out in PA, the "minimum" of "beyond 
the break" stars to the south would look similar to that of the 
north. Overall, the stars beyond the break radius appear to 
be spread out more along the major axis than the minor axis, 
and have the appearance of a nascent tidal tail system. It is 
difficult to understand how the observed increasing elliptic - 
ity of the system would originate and be sustained in a dSph 
with an extended DM halo; a very elongated luminous "halo" 
population with anisotropic velocity ellipsoid structure over- 
lapping and extending beyond a second, more concentrated 
population with a rounder distribution would be implied. 

(2) Our "King+break" radial density profile for Leo I is sim- 
ilar to the profiles observed in previous models of disrupting 
satellites (Johnston, Sigurdsson, & Hernquist 1999; Johnston, 
Choi, & Guhathakurta 2002, see also our new N-body sim- 
ulation presented in §6.2) as well as in the Sagittarius dSph 
(Majewski et al. 2003), an acknowledged example of a tidally 
disrupting satellite. 

(3) We have found a flat overall velocity dispersion profile 
for Leo I, a profile that is expected in the case of a tidally 
disrupting dSph (Kroupa 1997; Munoz et al., in preparation) 
and is actually observed in the case of the tidally disrupt- 
ing Sgr dSph (Majewski et al. 2004). While such profiles 

25 Careful examination of the stellar distribution in model galaxies for tidal 
disruption models show that stars released during the most recent perigalac- 
tic pass form a tidal tail-like feature mostly aligned with the semi-major axes 
of the satellite, whereas those that become unbound during older perigalac- 
tic passes tend to be spreaded out irrespective of the semi-major axes. One 
demonstration can be found in Figure 23, an orbital plot of our W-body sim- 
ulations where we use different colors for stars that become unbound dur- 
ing different perigalactic passes. Consequently, photometrically chosen dSph 
candidates located on degree scales away from the satellite's main body will 
not necessarily show a clearly defined tidal tail-like configuration. For this 
reason, we have limited our sample in r ma j when constructing the azimuthal 
sector count plot. 
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could also be formed in the presence of an extended DM halo 
(Kleyna et al. 2001), the detailed velocity distribution at large 
radius is not what is expected for the latter model, whereas it 
compellingly resembles (Fig. 19) a dynamically cold (likely 
bound) stellar population at the systemic Leo I velocity with 
a superposed population of unbound stars with larger velocity 
dispersion. 

(4) The asymmetry of the RVs of the outliers would seem 
to challenge the bound, extended halo population hypothesis, 
since in this model stars at a variety of orbital phases (i.e. 
both approaching and leaving the apocenters of their internal 
orbits) would be expected to yield a symmetrical velocity dis- 
tribution. On the other hand, such asymmetries are not only 
accommodated, but expected in tidal tail models (as we show 
in §6.2). 

We conclude that the most straightforward and natural in- 
terpretation for the observed Leo I structure and dynamics is 
that (1) it is tidally limited near the observed King limiting ra- 
dius with the dwindling bound population creating a dynami- 
cally cold signature at its outermost extent, and (2) increasing 
numbers of unbound stars are being observed at larger radius 
and forming nascent tidal tails that contribute the RV outliers 
at large radius. If this proposed model for the structure of 
Leo I is true, we may infer several things immediately about 
the mass loss history and orbit of this satellite galaxy, even 
without resorting to new N-body simulations. 

6.1.2. Inferred Mass Loss Rate 

Models of tidally disturbed satellites around the MW by 
Johnston, Sigurdsson, & Hernquist (1999) and Johnston, 
Choi, & Guhathakurta (2002) predict breaks in satellite ra- 
dial density profiles similar to that shown in Figure 9b. Under 
the assumption that Leo I stars found past the radial density 
profile break are unbound, extratidal debris, we may estimate 
the mass-loss rate. Johnston, Sigurdsson, & Hernquist (1999) 
provide an algorithm for determining the approximate mass- 
loss rate of a satellite using the strength of the break popula- 
tion. If the density profile of a satellite galaxy shows a break 
at rt,reak and the extra-break population is well defined out to 
a radius of r xt , the fractional mass-loss rate can be estimated 
from 
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where 9 is the angle between the line of sight and the plane 
perpendicular to the satellite velocity, g(9) is a geometric fac- 
tor corresponding to the orbit (which can be approximated 
as cos 0), « xt is the number of extratidal stars between 

^break 

and r xt , «break is the number of stars within rb rea k, and r ol -b is 
the period for a circular orbit at the present satellite Galacto- 
centric distance. For now we adopt the values of 9 and T or b 
for Leo I listed in Table 4 of Johnston, Sigurdsson, & Hern- 
quist (1999) (but provide a new model for this orbit in §6.2). 
By examining the radial density profile in Figure 9b, we de- 
termine rbreak = 10.2 arcmin. For r xt , we have adopted the 
major axis radial distance of the farthest data point in Fig- 
ure 9b since the corresponding surface density is near the 
background level. Applying corrections to account for miss- 
ing catalog giant stars falling in the Mosaic CCD chip gaps 
gives Hbreak = 1150. For the extratidal count, we scale our 
counts of stars at each annulus by the ratio of the total el- 
liptical area to the amount of that annulus in our actual field 
coverage, and we find n xt = 129. Using these values we obtain 
a mass-loss rate of (d//df)i = 3.1 x 10~ 3 Gyr" 1 , which implies 



an average of less than 1% of the total mass of Leo I is lost per 
gigayear. We note that the mass-loss rate equation (3) above is 
technically for satellites with extratidal radial density profiles 
that follow S xt ~ r~ l . In our case, the measured outer den- 
sity slopes are steeper than 7 = -1 (see Figure 9b). Johnston, 
Sigurdsson, & Hernquist (1999) also note that this equation is 
only good to within a factor of two. Therefore, the combined 
uncertainty in the mass-loss estimation is likely to be large. 
Nevertheless, the small net mass-loss rate for Leo I is consis- 
tent with the upper limit determined by Johnston, Sigurdsson, 
& Hernquist (1999), though our new density profile of Leo I 
has allowed us to re-calculate this improved upper limit to the 
mass-loss rate. 

For satellites with a less well-defined extratidal population, 
Johnston, Sigurdsson, & Hernquist (1999) derived an expres- 
sion for estimating the upper limit of the fractional mass- 
loss rate, and using more realistic models, Johnston, Choi, 
& Guhathakurta (2002) update the relation to: 

df\ ( 2 ~\ S xt (r bre ak) 7T _ 2 ,a\ 

-27rr bre ak , (4) 



dr 



^break 



where S xt (rb re ak) is the number density at r = rb rea k and r c i rc 
is same as T or \,. Using the values for the repeated parameters 
from above, we obtain an upper limit of (d//df)2 = 5.1 x 10~ 2 
Gyr" 1 , similar to what Johnston, Sigurdsson, & Hernquist 
(1999) found. 

Our results confirm the results of Johnston, Sigurdsson, & 
Hernquist ( 1 999) that among the dSphs studied by IH95 , Leo I 
apparently has one of the smallest mass-loss rates. The large 
orbit and few perigalactic passes of Leo I may explain this 
low mass-loss rate. For example, Leo II which has a lower 
RV at a similar distance to Leo I (which implies a smaller 
apogalacticon), has a higher estimated mass-loss rate. Despite 
the relatively low fractional mass loss rate derived for Leo I, 
that is has a perceptible one at all, and so clearly in the form 
of tidal tails, shows that no MW dSph with a smaller orbit and 
mass (i.e., virtually all other known Galactic dSphs) is likely 
to be immune from significant tidal effects (Johnston, Choi, 
& Guhathakurta 2002). 

6.2. A Tidally Disrupting Satellite Model to Explain the 
Outer Leo I Structure and Dynamics 

6.2.1. Fitting a Leo I Model 

Johnston, Choi, & Guhathakurta (2002) have A^-body sim- 
ulated a number of examples of disrupting satellite galaxies 
in MW-like potentials, including one system with mass and 
orbital properties likely similar to Leo I (their models 4 and 
5), and show how tidal disruption can be expected even in 
extreme cases of a satellite in a large orbit like that of Leo I. 
Here we undertake additional modeling specifically guided by 
the new Leo I observations to see whether we can obtain a 
somewhat closer model match to the newly observed Leo I 
structure and velocities found here. Our primary goal is to 
understand the nature of the extended population observed in 
the Leo I profile. More specifically, we seek to test whether 
a tidal disruption scenario can explain the observed (1) elon- 
gated dSph structure at radii comparable to and larger than the 
limiting radius (i.e. the bipolar distribution of the break popu- 
lation), (2) shape of the overall density profile, (3) dispersion 
profile over all radii, and (4) the asymmetric RV distribution. 
A disrupting model that simultaneously accounts for all or 
even some of these features would not only lend support to a 
tidal disruption scenario, but (as we shall show) also provide 
the first real constraints on the orbit of Leo I. 
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The iV-body simulation code we adopt has the same her- 
itage as the Johnston, Choi, & Guhathakurta (2002) mass- 
follows-light (i.e. single component) models. The static MW 
potential has a logarithmic, spherical (q = 1.0) halo with cir- 
cular speed 210 km s" 1 . The assumed solar distance to the 
Galactic center is 8.5 kpc and the total Galactic mass within 
50 kpc is 4.5 x 10 11 M Q . Other aspects of the potential are as 
in Law, Johnston, & Majewski (2005). 

The satellite is modeled by 10 5 particles originally config- 
ured as a Plummer (1911) model, 

_ GM sate ni tei o 



which has a physical scale length parameter, ro. The model is 
constrained so that the present satellite position and RV match 
those observed for Leo I: (l,b) = (226, +49)°, a 259 kpc helio- 
centric distance, and Galactocentric velocity (vgsr) of 180 km 
s -1 . However, because the proper motion of Leo I is unknown, 
this is a free parameter that ultimately determines the shape of 
the model Leo I orbit. After assuming a given present proper 
motion, the Leo I orbit (with the satellite as a point mass) is 
evolved back in time long enough to derive the phase space 
position of the satellite two apogalactica ago. The point mass 
is then replaced by the Plummer model satellite at this phase 
space position (after the satellite has been allowed to evolve 
and relax at infinite distance), and a full A^-body simulation is 
evolved forward to the present time. In general, for the mod- 
els tested here the start of the simulations occurs about 1 1-12 
Gyr ago. In the course of the modeling efforts for Leo I, some 
100 differently configured models have been run. 

Initial models were run to fix a likely orientation of the or- 
bit, under the assumption that the east-west position angle of 
the Leo I ellipticity and the orientation of the "break popu- 
lation" are caused by tidal effects and tidal stripping, respec- 
tively. jV-body simulations were run with the satellite hav- 
ing an orbital pole every 45° along its allowed set of poles 
(its "great circle pole family"; see Palma, Majewski, & John- 
ston 2002). It was found that a satellite with orbital pole near 
(l,b) = (122, +13)° yielded tails with the proper orientation, 
though this simulation cannot be discriminated (on the basis 
of the direction of Leo I's stretching alone) from one with a 
satellite having an opposite orbital pole (i.e. [302, -13]°). We 
adopt either of these antipodal orbital poles for the remain- 
ing simulations, which thereby fixes the position angle of the 
present Leo I proper motion to <~ 75° or <~ 255° in Galactic 
coordinates. 

The remaining simulation variables we explore are the ini- 
tial satellite mass and scale (ro) — which determine the size 
and density of the satellite and therefore the degree of satel- 
lite disruption — and the magnitude of the proper motion 
— which determines the properties (eccentricity, peri-and 
apogalacticon) of the orbit and also affects the degree of satel- 
lite disruption. Given the extreme distance and RV of Leo I, it 
seems probable that Leo I has a rather elliptical orbit (Taylor, 
Silk, & Babul 2004). Thus, our modeling efforts centered on 
orbits with perigalactica ranging from 10 kpc to 50 kpc (and 
consistent with the observed RV of Leo I); despite the large 
variation in perigalactica, this actually corresponds to orbital 
eccentricities ranging only from 0.80 to 0.96. With such orbits 
Leo I has an orbital period of about 6 Gyr. It is found that after 
fixing the satellite mass to of order that found in §4.2, a scale 
of ro = 0.3 ±0.1 kpc yields a final Leo I satellite with a tidal 
radius of order the observed Leo I King limiting radius. Thus, 



we adopt models with this general structure, and explore how 
varying the orbit shape is reflected in the resultant radial light 
profile and velocity distribution. For the latter, we "observe" 
radial velocities sampled from the model distribution in a spa- 
tial "footprint" mimicking that of the Keck DEIMOS spectra 
as shown in Figure 16; this gives the most direct and fair com- 
parison to the observed RV distribution seen in Figures 14 and 
18. 

It is found that a variety of simulations of Leo I on an ec- 
centric orbit can reproduce not only the overall "King+break" 
radial light profile that is characteristic of disrupting satel- 
lites (Johnston 1998), but also an asymmetric RV distribu- 
tion and a more or less flat velocity dispersion profile like that 
observed. To fine tune the model, we took account of three 
general correlations in turn (see Munoz, Majewski, & John- 
ston 2007): (1) The central velocity dispersion directly re- 
flects the adopted satellite mass. (2) The size of the observed 
"King profile" part of the satellite is set by the initial Plum- 
mer model scale ro. (3) With mass and scale set by the previ- 
ous conditions, the mass loss rate (hence the size of the break 
population) is then only driven by the orbital shape (i.e. im- 
pact parameter to the MW). Following these general guides, 
a narrow set of mass, scale and orbital shapes were found to 
give reasonable matches to the observed Leo I properties. The 
density and RV properties of two of these "matching" models 
are compared against the data in Figure 22, and their three- 
dimensional orbits are shown in Figure 23. 

The parameters and results for adopted models are listed in 
Table 5. Models 111 and 117 have been run in the MW poten- 
tial for 12.0 Gyr and 1 1.8 Gyr, respectively. In both cases the 
net mass loss is modest, with only 16% and 23% of the initial 
mass being unbound after 12 Gyr. We note that the implied 
mass loss rates for both models roughly match the estimates 
made in §6.1.2. As found in previously published simula- 
tions of satellite disruption (e.g., Law, Johnston, & Majewski 
2005), the break populations in the extreme orbit, Leo I model 
satellite density profiles are also found to be due to tidal dis- 
ruption; this is shown in Figure 22, where stars that are still 
bound as well as those that have been lost and become un- 
bound on each of the last two radial orbits are marked with 
separate symbols and colors. Thus, our original supposition 
that the observed Leo I break population may be due to un- 
bound stars is supported by the models. However, the relative 
density at which these breaks occur does not match that of the 
Leo I density profile if the beginning of the unbound debris 
is associated with the inner of the two radial profile breaks. 
Said another way, these models provide a good match to all 
the Leo I data if the first break in the radial density profile 
is considered as statistically insignificant or structurally unre- 
lated to unbound tidal debris. Even if the latter suppositions 
are not valid, we contend that these models still provide as 
good or a better simultaneous match to the overall observed 
structure and dynamics of any specific dSph than has been 
offered before. 

We also attempted to find a model that could accommo- 
date the inner density law break. Since the relative density 
at which tidally induced break populations occur is directly 
related to the instantaneous rate at which debris is being gen- 
erated, in order to match the relative position compared to the 
central density of the inner break at higher density we need 
to increase the mass loss rate of the satellite. This can be 
done by decreasing the satellite density. Figure 24 shows the 
properties of model 122 which has an initial mass of 6.5 x 10 7 
M©, a scale of 0.55 kpc and apo:perigalacticon of 450: 10 kpc. 
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Whereas the position of the model break is now matched to 
the inner of the two observed density breaks, the model satel- 
lite is found to be three times larger than the actual Leo I, and 
the velocity dispersion as well as the overall density profile 
are no longer a good match. We therefore find model 122 to 
provide a less satisfactory description of Leo I than models 
111 and 117. 

The final satellite masses for Models 111 and 117 (see Table 
5) are consistent with the Leo I mass derived from core fitting 
in §4.2 within the errors. The implied current total Leo I mass- 
to-light ratio from the models is M/L = 3 . 1 -4.5(M/L) . This 
lower M/L than found in §4.2 implies a relatively small DM 
content, and would bring Leo I in line with the M/L typical 
of dE galaxies and even globular clusters. It is important to 
recognize that we have obtained a good match to the Leo I RV 
asymmetry with both of the preferred models even though this 
observed property was not used as a model constraint. That 
such a result comes naturally lends further confidence to the 
plausibility of our models. The RV asymmetry in the models 
arises from the long extension of the trailing tidal arms to- 
wards the inner Galaxy (and the observer) seen in Figure 23 
and for which the innermost parts have a significant projec- 
tion along the line of sight in Models 111 and 117. However, 
the degree of that projection, and therefore the velocity spread 
of the asymmetry is obviously tied with the eccentricity of the 
orbit. This can be seen, for example, by the results (Fig. 25) of 
a series of models run with the same mass and scale, but more 
circular orbits than shown in Figure 23. As Figure 25 shows, 
less eccentric orbits have less asymmetry in their observed RV 
distribution. In addition, as might be expected, these models 
give poorer matches to the radial density and velocity disper- 
sion profiles. In this way, the observed RV distribution can 
apparently directly constrain the detailed shape of the Leo I 
orbit, and by direct comparison to the observed RV distribu- 
tion, we deduce that Leo I has an orbital eccentricity similar 
to that shown in Figure 23. We also deduce the general direc- 
tion of the Model 1 1 1/1 17 orbits (i.e. the general direction of 
the proper motion) to be correct (i.e. to the east), since a satel- 
lite with a similar eccentricity orbit but opposite orbital pole 
yields the opposite RV asymmetry for our Keck spectroscopy 
mask footprint (Fig. 16). This is demonstrated by the model 
results shown in the right-most panel in Figure 25. 

Finally, because we have shown by the models that stars in 
the break population are well described as nascent tidal tails, 
and by the analysis in §6.1.1 we have shown the Leo I break 
population lies predominantly east- west, we conclude that the 
orbital pole of (l,b) = (122, +13)° is approximately correct. 
Thus, based on the apparent discriminatory power of the Fig- 
ure 22 parameters as well as the direction of the tidal arms, 
we conclude that the orbit shown in Figure 23 to be a reason- 
able approximation to the true Leo I orbit. A check on the 
orbit will obviously be delivered by the measurement of the 
Leo I proper motion, which can be expected in a decade or 
so from a key project (led by SRM) of NASA's Space Inter- 
ferometry Mission. The expected current proper motions for 
Leo I are predicted to be (/Lt/cosfr, Hb) = (0.0046,0.0219) and 
(0.0020,0.0138) mas yr" 1 , respectively for Models 111 and 
117 (including solar motion, which is assumed to have a 232 
km s" 1 revolutionary speed about the MW, and additional pe- 
culiar motion of 9.0 and 7.0 km s" 1 in the Galactic radial and 
Z directions, respectively). 

In the meantime, verification of the overall picture painted 
by our model fitting here, including the orbit, can be made by 
obtaining more RVs on the east side of Leo I. Our models pre- 



dict that at large radii on the east side of Leo I an asymmetry of 
the RV distribution should be observed opposite (i.e., toward 
lower RVs) that we have found to the west. The change in the 
sense of asymmetry arises from the sampling of the leading 
arm of the Leo I tidal tails on the east side, whereas our west- 
ern RVs have been sampling trailing arm debris. Mapping the 
Leo I tidal tails to larger angular separations would also pro- 
vide significant leverage on the Leo I orbit and mass loss his- 
tory. Figure 23 suggests that quite long tails should exist for 
two orbits of mass loss. However, the problem is quite obser- 
vationally challenging because the tails are at fairly low sur- 
face brightness and the corresponding densities of Leo I tidal 
tail giant stars will eventually become quite sparse. A search 
for the main sequence turn-off CMD feature for the Leo I tidal 
arms near the satellite will require reaching to V > 28 over 
large areas. From the ground, such work is severely hampered 
by the difficulty of star/galaxy separation at these magnitudes 
(D. Martinez-Delgado, private communication). 

6.2.2. Implications for the Local Group Path of Leo I and 
Mass of the Galaxy 

The provenance of Leo I has important implications for the 
mass of the MW. The high RV of Leo I at 259 kpc translates 
to a large implied lower limit to the escape velocity of the 
MW if this satellite is bound. If Leo I has made two orbits 
about the MW, as assumed in our models, then it must be 
bound. If it has made only one perigalactic passage, Leo I 
could still be bound to the MW, having become so on the last 
orbit; however, in this case Leo I may also be unbound and on 
a hyperbolic orbit. 

Unfortunately, it is not clear that we can, with certainty, 
establish whether two or only one orbits have occured. Ac- 
cording to the models, much of the unbound debris observed 
near Leo I has detached in the last orbit (see Fig. 23). To ver- 
ify whether any differences can be discerned between a one 
and two orbit Leo I, we ran a simulation (Model 118) with the 
same satellite and orbital properties as Model 111, but where 
only one radial orbit has occured (i.e. the model is started 4.3 
Gyr ago). Figure 26 shows that the radial density profile from 
such a simulation is virtually indistinguishable from the two 
orbit Model 111. This is because of the small contribution 
of older debris to the density of Leo I over the currently ob- 
served area. On the other hand, this small contribution does 
become more obvious in its influence on the observed veloc- 
ity dispersion, as seen in the bottom panels of Figure 26: The 
one orbit model has a smaller velocity dispersion at large radii 
compared to that of the two orbit model, where older debris 
helps inflate the dispersion. The latter, two orbit dispersion 
profile is a closer match to that observed for Leo I, and so, 
based on this evidence alone can we tentatively suggest that a 
two orbit model is favored for Leo I. 

The star formation history of Leo I may lend circumstantial 
support to the two orbit scenario. Models 111 and 117 sug- 
gest that Leo I endured fairly substantial tidal shocking both 
<~ 1 and ^7.5 Gyr ago. This general orbital picture seems 
to coincide with the star formation history of Leo I studied by 
Gallart et al. (1999b), who finds that most of its star formation 
activity occurred between 7 and 1 Gyr ago. The oldest (> 10 
Gyr) Leo I population, discovered by Held et al. (2000, blue 
HB stars) and Held et al. (2001, RR Lyrae variables), likely 
formed as a result of the initial collapse of gas that led to the 
formation of Leo I. The "major" star formation that started 
~ 7 Gyr ago is roughly timed with the previous perigalactic 
passage in our model (see also models by Mayer et al. 2001). 
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We can only speculate as to the causes of the abrupt drop in 
star formation at <~ 1 Gyr ago. Leo I may have simply ran 
out of gas to fuel the star formation. The absence of gas in 
Leo I 26 (Knapp, Kerr, & Bowers 1978; Bowen et al. 1997) is 
consistent with this picture. The coincidence that Leo I passed 
through the inner Galaxy about this time (and possibly with 
a rather small impact parameter; Fig. 23) may also suggest 
gas stripping by the MW disk or perhaps by a dense high ve- 
locity cloud. A massive (3 x 10 6 M Q ) fragmented H I cloud 
structure around Leo I (Blitz & Robishaw 2000) perhaps indi- 
cates that the gas has been perturbed by tidal shocking. Nev- 
ertheless, something triggered star formation at a time that 
matches well the timing of the first perigalacticon in our two 
orbit model, and, together with the velocity dispersion data, 
our results support the hypothesis of Z89 that Leo I is bound 
to the MW and has had two radial orbits in this state. 

If Leo I is bound, a large MW mass is implied. The total 
mass of the MW within 260 kpc from the center in our simu- 
lations is 1.8 x 10 12 Mq, and it has a profile yielding a mass 
interior to 50 kpc that is consistent with those found in the 
analyses by Z89, Wilkinson & Evans (1999), and Sakamoto, 
Chiba, & Beers (2003) when these analyses include Leo I as 
a bound satellite. Verifying the length of the Leo I tidal tails 
would verify whether it has had multiple orbits around the 
MW and is thus bound. Obviously, a definitively measured 
proper motion for Leo I would also help determine the true 
orbit of Leo I and refine the mass of our Galaxy. 

What do our results say about the origin of Leo I? First, it is 
probable that a hyperbolic orbit for Leo I would produce sim- 
ilar results (over the sky area we have surveyed here) to those 
seen for a single orbit, bound model. Moreover, such an orbit 
will have an overall shape not dissimilar in overall orientation 
and general direction to the last orbit of Models 111 and 117, 
and thus we may determine from which direction Leo I ap- 
proached the MW in this case. Interestingly, in either the one 
orbit (bound or hyperbolic) or two orbit case, our modeling 
seems to rule out any close association of Leo I with M3 1 over 
relevant timescales to tidal stripping: (1) In the bound case, 
Figure 23 shows that Leo I is not in the vicinity of M3 1 (which 
has current MW coordinates of (X , Y, Z) cc ~ (375 , 620,-285) 
kpc) since approximately two orbital periods ago, and even 
then the inferred distance between M31 and Leo I is <~ 700 
kpc (and this ignores the motion of M3 1 over a Hubble time!). 
(2) In the hyperbolic case, Leo I would have apparently en- 
tered the Local Group in an orbital plane almost perpendic- 
ular to the direction of M31 (i.e., from the general direction 
of the last apogalactica shown in Figure 23). Were the Byrd 
et al. (1994) hypothesis that Leo I was once bound to M31 to 
be true, it would have had to have been released from M3 1 at 
least ~ 10 Gyr ago. 

Finally, we add a note of caution about interpreting the po- 
tential shape of the Leo I orbit as inferred or extrapolated from 
Figure 23, which represents a model run in a static MW po- 
tential. Obviously, currently favored hierarchical galaxy for- 
mation models imply a continuing growth of the Galactic po- 
tential with time, and this will alter the orbits of satellites. 
On the other hand, the MW formation models of Bullock & 
Johnston (2005) suggest that destroyed satellites contributing 
mass to the MW halo accrete predominantly before the last 
two postulated Leo I perigalactica (i.e., before 7 Gyr ago). 
Thus, at least over this timescale, a static MW approximation, 

26 (Blitz & Robishaw 2000) finds extended H I emission around Leo I but 
not toward the galaxy itself. 



and the resulting Figure 23 orbit, may be appropriate descrip- 
tors. Moreover, as recently shown by Penarrubia et al. (2006), 
tidal streams respond adiabatically to evolving potentials, so 
that even if the MW potential evolved over the past 7 Gyr or 
so, the Leo I tidal stream visible today would resemble that 
derived in a static MW potential with the present day mass 
profile. 

7. SUMMARY AND CONCLUSIONS 

We have photometrically surveyed a 4.5 degree 2 region cen- 
tered on Leo I in the M, Ti, DD051 filter system in order 
to explore the extended morphology of this dwarf spheroidal 
galaxy, which is currently regarded as the most distant of the 
known Galactic satellites (unless the more distant Phoenix 
system is bound to the MW). 

The photometric data were used to select Leo I giant can- 
didates based on the two criteria: (1) the gravity sensitive 
(M-DD051, M- Tj) diagram, which separates distant giant 
stars from contaminating, foreground, metal-rich disk dwarfs, 
and (2) CMD positions commensurate with the temperature- 
apparent magnitude combination of stars on the Leo I giant 
branch. The background level of our "Leo I giant star" sample 
is determined to be small, and a 100% reliability in the iden- 
tification of bona fide Leo I giants is found via testing with a 
total of 133 stars in the Leo I field with previously published 
or new Keck spectroscopy. 

We derive a new set of Leo I structural parameters by fitting 
a single-component King profile to the Leo I giant candidates. 
Coupling this to the central velocity dispersion we have mea- 
sured (§5.5), we use core-fitting techniques to derive a mass 
for the Leo I system of 3.5 x 10 7 M Q and a total (M/L) v of 
4.6M Q /L y, 27 values not too dissimilar from previous study 
using the same technique (M98). The two-dimensional dis- 
tribution of Leo I giant candidates shows many giant stars 
outside the derived King limiting radius. These are primarily 
along the major axis and spectroscopy of a subsample of these 
"extratidal" stars shows they are actually associated with the 
dSph. This population of stars shows up as a break at a ma- 
jor axis radial distance of <~ 10 arcmin in the radial density 
profile. Our new Keck spectroscopy confirms this excess of 
stars beyond the King limiting radius to be made of true Leo I 
members. 

Our Leo I velocity dispersion profile shows a flat and then 
rising trend at large radii. We also find that the Leo I RV 
distribution, particularly for stars at large angular separations 
(which in our data are primarily to the west side of Leo I), in- 
cludes a population of stars with a broad and skewed, asym- 
metrical distribution toward positive RVs overlapping a sec- 
ond, very much colder population at the nominal Leo I mean 
velocity. 

We interpret these features as support for a picture wherein 
Leo I has been tidally disrupted on at least one, but at most 
two, perigalactic passages by a massive Local Group mem- 
ber (most likely the MW for both), and find these phenomena 
naturally produced in mass-follows-light TV-body simulations 
of a tidally disrupting Leo I analog in a Milky Way-like po- 
tential. The best-matching TV-body simulations to both the 
observed structural morphology and velocity distribution of 
Leo I are those where the tidally disrupting satellite is on a 
rather high eccentricity (0.93-0.96), small perigalacticon (10- 
15 kpc), bound orbit around the Milky Way, and has a present 

27 These values are 14% smaller than those in §4.2, where the M98 central 
velocity dispersion was used to derive the mass. 
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total M/L <~ 3 -4.5(M/L) . These best fitting model masses 
are 58-85% the M/L derived from the central velocity disper- 
sion and core fitting, but it is not unreasonable to presume 
that the latter method yields somewhat inflated masses be- 
cause of systematic increases in the true dispersion due to as- 
trophysical processes such as the presence of binaries and the 
atmospheric jitter common to giant star atmospheres. Thus, 
the likely M/L of the satellite is rather modest and not un- 
like those of other elliptical systems of approximately similar 
mass scale that are typically regarded as low in, or devoid of, 
DM (dE galaxies and globular clusters). 

Because of the rather close match between our model re- 
sults and observations, and because disrupting satellites on 
highly radial orbits appear to yield great discriminatory power 
in this regard, we have been able to constrain the likely orbit 
of Leo I without the measurement of its proper motion. The 
orientation of the current satellite orbital plane can be fixed by 
matching model tidal tails to the predominant direction of the 
observed break population (i.e., more or less along the major 
axis), whereas we find that the direction of angular momen- 
tum in this orbital plane is well constrained by the sense of the 
RV asymmetry we have observed to the west side of the satel- 
lite. Our models demonstrate that a positive RV asymmetry in 
the models is produced by trailing tidal debris for the reced- 
ing satellite. Thus, we predict that an opposite RV asymmetry 
will be found on the east side of Leo I from leading tidal de- 
bris. Such a result would provide an important verification of 
the tidal disruption model we have put forward in this paper. 

Our observed Leo I RV distribution is most consistent with 
a two MW orbit history for Leo I, with both orbits around the 
MW, however we cannot yet definitively rule out a one orbit 
scenario. However, whether Leo I is bound to the MW or 
on an unbound, hyperbolic orbit around the MW, our results 
seem to rule out a Leo I orbit that includes a previous associ- 
ation with M31 within the last 10 Gyr, in contradistinction to 
the Byrd et al. (1994) scenario of a relatively recent origin of 
Leo I from M3 1 . 

Leo I has long played a "spoiler" role in setting the mass of 
the MW because of its huge rv 2 lever arm in Jean's equation- 
based determinations using tracer particles of unknown proper 
motion (i.e. unknown orbit size and shape. The large im- 
plied escape velocity at Leo I's distance implies large MW 
masses if Leo I is bound to our galaxy (Z89, Wilkinson & 
Evans 1999), although its influence on the mass determina- 
tion is lessened when more complete samples of objects with 
complete phase space data are employed (e.g., Sakamoto, 
Chiba, & Beers 2003). Our observed RV distribution for Leo I 
slightly favors a two perigalacticon pass, bound MW orbit for 



Leo I, thereby suggesting that higher mass estimates for the 
MW may be more correct. 

In general, our tidally disrupting mass-follows-light satel- 
lite models provide a quite satisfactory match to the observed 
properties of Leo I, but a few details — namely the apparent 
double break in the density profile and the density at which 
the inner of these breaks occur — we have yet to account for 
in these initial modeling efforts. Nevertheless, we contend 
that the scenario of a tidally disrupting, low M/L system on 
a highly radial orbit provides a rather complete explanation 
for the observed properties of Leo I. While some properties 
of Leo I and other dSphs (e.g., flat velocity dispersion profiles 
and break populations) have also been explained by postulat- 
ing that these systems are embedded in extended DM halos, 
such an explanation in the case of Leo I appears less com- 
pelling in that it cannot account for as many of the observed 
properties of the system (e.g., the asymmetry in the RV dis- 
tribution and increasing ellipticity with radius). In contrast, 
tidal tails by now provide a well-established observational 
paradigm for dSph satellites, with the Sagittarius dSph the 
most vivid example. We contend that the tidal disruption of 
Sagittarius is not unique and that Leo I may be another ex- 
ample of the phenomenon, albeit at a much lower mass loss 
rate commensurate with a satellite on such an extreme orbit. 
Tidal disruption observed among both the closest and farthest 
of the Milky Way satellites suggests that this process may be 
ubiquitous, and that similarly structured satellites with ellip- 
tical orbits and distances between these two examples might 
also be expected to be experiencing tidally induced mass loss, 
likely with an intermediary range of mass loss rates. 
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TABLE 1 
LEO I OBSERVATION LOG. 





RA 


Dec 


Exp time( 


;ec) / Airmass / FWHM(arcsec) 


Field 


(J2000) 


(J2000) 


M 


1 2 


DD051 




1 0-08-28 69 


+12-18-17 2 


70/1 16/12 


70 / 1.15 / 0.8 


700 / 1.13 / 0.9 


NN 


1 0-08-29 27 


+13:18:16.8 


70 / 1.20 / 1.2 


70 / 1.19 / 0.9 


700 / 1.17 / 1.0 


EE 


10:12:28.46 


+12:18:08.7 


70/ 1.29 / 1.1 


70 / 1.28 / 0.9 


700 / 1.25 / 1.0 


WW 


10:04:28.91 


+12:18:25.2 


70 / 1.46 / 1.2 


70 / 1.44 / 1.0 


700 / 1.39 / 1.1 


N 


10:08:28.12 


+12:48:23.4 


100 / 1.08 / 2.2 


200 / 1.09 / 1.8 


1000 / 1.10 / 2.3 


E 


10:10:28.12 


+12:18:23.4 


100 / 1.21 / 2.0 


200 / 1.23 / 1.6 


1000 / 1.27 / 1.8 


W 


10:06:28.11 


+12:18:23.4 


100 / 1.15 / 1.8 


200 / 1.16 / 2.0 


1000 / 1.19 / 1.7 


NE 


10:10:28.12 


+12:48:23.4 


100 / 1.35 / 1.6 


200 / 1.38 / 1.6 


1000 / 1.44 / 1.7 


NW 


10:06:28.55 


+12:47:57.1 


100 / 1.13 / 1.9 


200 / 1.14 / 2.1 


1000 / 1.20 / 1.7 


EEE 


10:14:28.12 


+12:18:23.4 


100 / 1.09 / 1.1 


200 / 1.09 / 1.2 


1000 / 1.11 / 1.3 


NEE 


10:12:28.11 


+12:48:23.4 


100 / 1.13 / 1.4 


200 / 1.14 / 1.4 


1000 / 1.16 / 1.3 


NEEE 


10:14:28.12 


+12:48:23.4 


100 / 1.19 / 1.3 


200 / 1.20 / 1.5 


1000 / 1.24 / 1.8 


NEEEE 


10:16:28.12 


+12:48:23.4 


100 / 1.27 / 2.5 


200 / 1.29 / 2.2 


1000 / 1.34 / 1.8 


WWW 


10:02:30.17 


+12:18:23.4 


100 / 1.76 / 1.5 


200 / 1.81 / 1.5 


1000 / 1.95 / 1.4 


SWW 


10:04:28.12 


+11:48:23.4 


100 / 1.26 / 1.6 


200 / 1.27 / 1.9 


1000 / 1.32 / 1.6 


SWWW 


10:02:28.12 


+11:48:23.4 


100 / 1.48 / 1.8 


200 / 1.51 / 1.8 


1000 / 1.59 / 1.4 



TABLE 2 
Number of 
color-color selected 

giants in rgb 
bounding boxes offset 
along brighter 
magnitudes. 



AM 


Number of Giants 


0.0 


1196 a 


-0.5 


709 a 


-1.0 


253 a 


-1.5 


27 a 


-2.0 


5 


-2.5 


2 


-3.0 


1 


-3.5 


5 


-4.0 


4 


-4.5 


4 


-5.0 


b 


-5.5 


o b 



"Sample contains stars 
from the Leo I system. 

b Sample has an incomplete 
sampling of the CMD at the 
bright end. 



TABLE 3 
Fit Parameters for Leo I. 







r, 


e 


P.A. a 


Reference 


(arcmin) 


(arcmin) 


(=1 -bid) 


(deg) u b 


This study (King) 


5.4 ±0.4 


13.4±0.7 


0.37 ±0.02 


84±2 


This study (PLC) 


7.3 ±0.8 




0.37 ±0.02 


85 ±2 4.8 ±0.6 


Hodge 1963 




14.3 ± 1.0 


0.31 ±0.07 




Hodge 1971 


4.5 


13.9±0.5 


0.31 ±0.07 




IH95 


3.3 ±0.3 


12.6±1.5 


0.21 ±0.03 


79±3 



"Major-axis position angle measuring from North=0° to East=90° . 
b See equation (1) of Kleyna et al. (1998). 



22 



TABLE 4 

Observed Radial Velocities of Stars 



Slit 




RA 


Dec 


Vhel* V GS R b 


Tell. Off. c 


No. 


ID 


(J2000.0) 


(J2000.0) 


M (M-T 2 )o (M-D51) (kms" 1 ) (kms" 1 ) 


(kms- 1 ) CCP d g e 



KD 1-00 


A1 7R9 f 

Ol / OZ 


1U.U / .jj.jj 


-4-17-7fv70. 1 

Tli.iO.iU. 1 


1 8 6.Q 

lo.DV 


1 

1 .OO 


—0 7^ 


U.J 


1UU.O 


—41 4 


VI. Jy 


j 


K U 1 -U 1 


ol /3o 


10.0/. 40. 39 


i 1 7.7A-T7 H 
+ 12.20.2 1 . 1 


io n 
lo.l / 


1 

1.65 


n on 
—0.20 


39.3 


— o/.o 


o c n 
35 .9 


a Art 
0.40 


6 


KD1-02 


72520 f 


10:07:30.84 


+12:13:41.8 


17.99 


0.79 


0.03 


158.6 


50.7 


— 17.8 


0.39 


4 


KD1-03 


72418' 


10:07:44.50 


+12:12:36.8 


16.44 


0.81 


0.04 


-13.1 


-121.0 


-18.9 


0.50 


5 


KD1-04 


61774 


10:07:34.30 


+12:25:09.4 


18.93 


2.90 


-0.04 


62.5 


-44.7 


-16.0 


0.46 


5 


KD1-05 


61788 


10:07:31.61 


+12:23:30.4 


17.56 


1.46 


-0.17 


26.2 


-81.1 


-23.5 


0.59 


5 


KD1-06 


61770 


10:07:34.86 


+12:25:15.0 


18.77 


2.34 


-0.17 


12.4 


-94.8 


-22.8 


0.57 


6 


KD1-07 


61785 


10:07:32.43 


+12:24:50.0 


19.84 


3.14 


-0.04 


10.9 


-96.3 


-30.7 


0.34 


4 


KD1-08 


72466 


10:07:38.44 


+12:13:33.4 


20.00 


2.84 


-0.02 


24.0 


-83.9 


-14.5 


0.53 


6 


KD1-09 


61692 


10:07:45.33 


+12:23:47.6 


20.56 


2.96 


-0.08 


4.0 


-103.3 


-23.9 


0.50 


6 


KD1-10 


72449 


10:07:40.94 


+12:13:03.6 


20.11 


1.82 


0.02 


288.3 


180.4 


-8.0 


0.91 


6 


KD1-11 


72472 


10:07:37.98 


+12:17:44.3 


20.19 


1.65 


0.01 


293.5 


185.9 


-14.0 


0.87 


6 


KD1-12 


72456 


10:07:39.52 


+12:16:54.6 


20.24 


1.58 


0.04 


303.1 


195.4 


-9.5 


0.64 


7 


KD1-13 


61733 


10:07:40.43 


+12:20:23.4 


20.19 


1.50 


0.03 


288.1 


180.6 


-11.2 


0.80 


6 


KD1-14 


72526 


10:07:30.06 


+12:16:21.6 


20.41 


1.48 


0.00 


290.1 


182.4 


-10.8 


0.82 


6 


KD1-15 


72493 


10:07:34.86 


+12:12:43.9 


20.57 


1.51 


0.00 


291.8 


183.9 


-10.4 


0.75 


6 


KD1-16 


72463 


10:07:38.60 


+12:11:28.7 


20.55 


1.47 


0.04 


278.9 


170.9 


-9.0 


0.74 


6 


KD1-17 


61790 


10:07:30.46 


+12:21:32.6 


20.65 


1.53 


0.02 


302.0 


194.6 


-12.9 


0.53 


6 


KD1-18 


61757 


10:07:37.01 


+12:22:42.9 


20.77 


1.48 


0.01 


289.1 


181.8 


-13.4 


0.65 


7 


KD1-20 


72407 


10:07:45.47 


+12:15:32.7 


21.31 


1.30 


0.02 


274.8 


167.1 


-5.7 


0.49 


7 


KD1-21 


72432 


10:07:43.00 


+12:16:27.9 


21.37 


1.33 


0.05 


293.4 


185.7 


-3.7 


0.49 


6 


KD1-22 


61783 


10:07:32.88 


+12:21:01.9 


21.42 


1.35 


0.01 


294.6 


187.1 


-14.6 


0.52 


6 


KD1-23 


61755 


10:07:37.37 


+12:22:27.3 


21.41 


1.32 


0.07 


313.9 


206.5 


-16.2 


0.36 


6 


KD1-24 


61753 


10:07:37.44 


+12:19:18.2 


21.49 


1.59 


0.22 


283.3 


175.8 


+3.8 


0.35 


6 


KD1-25 


61746 


10:07:38.47 


+12:19:27.3 


21.93 


1.26 


0.11 


288.1 


180.6 


-11.3 


0.28 


4 


KD1-26 


72464 


10:07:38.53 


+12:13:11.0 


21.56 


1.21 


0.06 


290.9 


183.0 


-9.0 


0.34 


6 


KD1-27 


72525 


10:07:30.36 


+12:11:44.5 


21.62 


1.27 


0.06 


314.8 


206.8 


-8.2 


0.25 


5 


KD1-28 


61769 


10:07:34.97 


+12:24:23.2 


21.40 


3.42 


0.07 


53.2 


—54.1 


-14.5 


0.42 


6 


KD1-29 


72505 


10:07:33.39 


+12:11:00.4 


18.90 


0.85 


0.03 


178.3 


70.3 


-9.7 


0.42 


5 


KD1-30 


61743 


10:07:39.55 


+12:26:39.4 


18.98 


0.81 


0.03 


129.7 


22.6 


-14.8 


0.28 


4 


KD1-32 


72490 


10:07:35.03 


+12:15:48.2 


19.88 


1.45 


-0.14 


52.7 


-55.0 


-6.3 


0.57 


6 


KD1-33 


61707 


10:07:43.58 


+12:21:21.2 


19.41 


0.91 


0.06 


219.1 


111.7 


-14.7 


0.18 


4 


KD1-34 


61694 


10:07:44.99 


+12:19:52.8 


21.70 


3.14 


0.03 


24.2 


-83.3 


-9.1 


0.47 


6 


KD1-36 


61766 


10:07:35.12 


+12:22:05.2 


20.88 


2.10 


-0.22 


68.6 


-38.8 


-5.9 


0.57 


6 


KD1-37 


61772 


10:07:34.39 


+12:18:57.0 


20.78 


1.76 


-0.24 


74.9 


-32.7 


+4.1 


0.48 


6 


KD1-38 


61775 


10:07:34.18 


+12:22:52.2 


20.72 


1.29 


-0.06 


-52.6 


-159.9 


-21.0 


0.48 


4 


KD1-39 


61724 


10:07:41.49 


+12:23:16.0 


21.44 


1.85 


-0.22 


24.9 


-82.4 


-19.2 


0.44 


5 


KD1-40 


72531 


10:07:29.23 


+12:11:59.5 


20.38 


0.76 


0.05 


318.1 


210.1 


-10.3 


0.18 


4 


KD1-42 


72439 


10:07:42.05 


+12:11:20.8 


21.21 


1.27 


0.00 


307.5 


199.5 


-13.7 


0.35 


5 


KD1-45 


72508 


10:07:32.88 


+12:15:58.1 


21.43 


1.22 


-0.01 


282.7 


175.0 


-1.1 


0.59 


7 


KD1-46 


72455 


10:07:39.82 


+12:17:55.1 


21.53 


1.31 


-0.04 


316.1 


208.5 


+6.7 


0.31 


4 


KD1-50 


61792 


10:07:30.07 


+12:22:20.6 


21.88 


1.13 


0.03 


305.3 


197.9 


-11.5 


0.34 


6 


KD2-03 


21673 


10:08:36.93 


+12:20:11.5 


19.47 


1.99 


-0.02 


299.0 


191.7 


-2.8 


1.08 


7 


KD2-04 


20868 


10:08:44.71 


+12:20:20.9 


19.79 


1.72 


0.04 


292.2 


185.0 


-1.9 


0.93 


6 


KD2-05 


23429 


10:08:25.15 


+12:20:07.5 


19.85 


1.78 


-0.02 


294.0 


186.7 


-4.9 


1.00 


6 


KD2-06 


31007 


10:08:40.95 


+12:18:06.7 


19.94 


1.83 


0.03 


293.3 


185.9 


-1.9 


0.88 


7 


KD2-07 


20898 


10:08:44.30 


+12:19:56.3 


19.89 


1.78 


0.03 


285.3 


178.0 


-4.1 


0.93 


7 


KD2-08 


20945 


10:08:43.74 


+12:19:57.9 


19.84 


1.70 


0.03 


296.9 


189.6 


-3.8 


0.98 


6 


KD2-09 


23059 


10:08:27.55 


+12:20:20.9 


19.81 


1.65 


-0.02 


276.9 


169.6 


-1.5 


0.80 


7 


KD2-11 


20198 


10:09:07.46 


+12:18:38.7 


19.99 


1.77 


0.01 


291.2 


183.9 


-5.0 


0.88 


6 


KD2-12 


20231 


10:09:05.15 


+12:20:14.4 


20.05 


1.73 


0.01 


288.8 


181.6 


-6.9 


0.94 


6 


KD2-13 


20345 


10:08:58.71 


+12:18:36.4 


20.07 


1.70 


0.00 


284.5 


177.2 


-4.8 


0.88 


7 


KD2-14 


20982 


10:08:43.32 


+12:21:23.6 


20.10 


1.72 


0.05 


286.0 


178.8 


-5.3 


0.80 


6 


KD2-15 


22396 


10:08:31.81 


+12:19:56.0 


20.00 


1.58 


-0.02 


287.5 


180.2 


+0.6 


0.70 


6 


KD2-16 


20129 


10:09:14.89 


+12:18:35.0 


20.10 


1.65 


0.01 


275.5 


168.3 


-8.3 


0.74 


7 


KD2-17 


20248 


10:09:04.17 


+12:20:46.0 


20.13 


1.68 


0.03 


292.0 


184.8 


-11.7 


0.93 


6 


KD2-18 


23131 


10:08:27.15 


+12:20:06.3 


20.09 


1.63 


0.05 


280.1 


172.8 


-5.4 


0.86 


6 


KD2-19 


31135 


10:08:39.39 


+12:17:23.8 


20.13 


1.67 


0.09 


286.7 


179.3 


+1.2 


0.81 


6 


KD2-20 


70295 


10:08:22.67 


+12:17:21.3 


20.21 


1.75 


0.00 


287.2 


179.7 


-6.4 


0.42 


5 


KD2-22 


21617 


10:08:37.29 


+12:20:12.0 


20.20 


1.68 


0.01 


291.9 


184.6 


-1.5 


0.90 


6 


KD2-23 


60019 


10:08:24.25 


+12:18:53.5 


20.20 


1.65 


0.03 


284.0 


176.6 


+4.2 


0.84 


6 


KD2-24 


30123 


10:09:13.77 


+12:17:23.0 


20.28 


1.71 


0.04 


296.9 


189.6 


-6.1 


0.75 


6 


KD2-25 


30522 


10:08:49.87 


+12:18:01.2 


20.21 


1.61 


0.04 


291.4 


184.0 


+0.9 


0.79 


6 


KD2-26 


60310 


10:08:21.89 


+12:19:33.4 


20.25 


1.63 


0.02 


280.4 


173.0 


+2.4 


0.79 


6 


KD2-27 


30102 


10:09:15.88 


+12:17:18.5 


20.39 


1.74 


0.03 


284.1 


176.8 


-3.8 


0.80 


6 


KD2-28 


21961 


10:08:34.69 


+12:19:11.4 


20.27 


1.56 


0.03 


283.6 


176.3 


-1.1 


0.82 


6 


KD2-29 


21005 


10:08:42.96 


+12:19:55.7 


20.30 


1.58 


0.02 


272.2 


164.9 


-3.6 


0.79 


6 


KD2-30 


32615 


10:08:27.92 


+12:17:44.7 


20.23 


1.45 


0.06 


279.6 


172.1 


+2.1 


0.82 


6 


KD2-31 


30710 


10:08:45.11 


+12:17:33.6 


20.45 


1.63 


0.02 


279.4 


172.0 


-0.9 


0.68 


6 


KD2-32 


20514 


10:08:52.57 


+12:19:05.9 


20.45 


1.58 


0.06 


302.5 


195.2 


-4.9 


0.78 


6 


KD2-34 


23305 


10:08:25.92 


+12:20:57.8 


20.33 


1.41 


0.03 


280.2 


172.9 


-6.7 


0.81 


6 


KD2-35 


23364 


10:08:25.55 


+12:18:46.6 


20.36 


1.44 


0.05 


276.8 


169.4 


-3.5 


0.72 


6 


KD2-36 


30941 


10:08:41.79 


+12:17:15.3 


20.44 


1.48 


0.03 


284.2 


176.8 


-4.2 


0.69 


7 


KD2-37 


60143 


10:08:23.12 


+12:18:55.3 


20.42 


1.45 


0.02 


301.3 


193.9 


+7.3 


0.65 


6 


KD2-39 


32487 


10:08:28.73 


+12:17:39.1 


20.38 


1.40 


0.09 


297.1 


189.6 


-2.3 


0.73 


6 


KD2-40 


30156 


10:09:10.01 


+12:18:05.2 


20.46 


1.46 


0.08 


291.1 


183.8 


-5.5 


0.61 


6 


KD2-41 


20491 


10:08:53.21 


+12:20:40.6 


20.54 


1.51 


0.04 


279.4 


172.2 


-4.9 


0.81 


7 


KD2-42 


22143 


10:08:33.48 


+12:20:34.5 


20.53 


1.49 


0.03 


288.8 


181.5 


-4.0 


0.82 


6 
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TABLE 5 

Parameters and results for adopted models. 



Model 


(10 7 M Q ) 


(kpc) 


'pen 
(kpc) 


r d 

i. apo 

(kpc) 


pe 

(Gyr) 


T f 
(Gyr) 


(10~ 2 Gyr 1 ) 


(1O 7 M ) 


111 
117 


3.0 
4.0 


0.3 
0.3 


15 
10 


450 
450 


6.33 
6.05 


0.96 
0.95 


1.3 
2.2 


2.3 
3.4 



a Initial mass of model satellite. 

k Scale length parameter in equation (5). 

c Perigalactic distance. 

^Apogalactic distance. 

e Radial period. 

^Time since last closest approach. 
^Fractional mass-loss rate. 
^Final mass of model satellite. 
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FIG. 2. — Maps of (a) all stars detected in our survey and (b) stars brighter than M = 21.5. The solid lines give rough indication of the boundaries of each 
observed field. The ellipses represent the Leo I tidal boundary derived by IH95 (r = i2'.6,PA = 79°; e = 0.21) and the asterisk symbols mark the positions of the 
1st magnitude star Regulus. The inhomogeneities in density among the fields of (a) is a reflection of variation in limiting magnitude across our survey area due 
to different observational conditions. Note that the inhomogeneities are absent in (b). 
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FIG. 4. — (M — T2, M — D5l)o diagrams for stars in (a) C frame and (b) all other frames. The box drawn with solid lines shows the bounding region we have 
employed to select metal-poor giant star candidates. 
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FIG. 5. — (M —T2,M)o color-magnitude diagram for stars selected as metal-poor giants in Figure 4 (a) within the tidal boundary derived by IH95 and (b) within 
our entire survey area. The bounding boxes shown by the solid lines are our CMD selection criteria. 
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FIG. 6. — Positions of the previously identified Leo I giant stars shown along with our catalog of stars with photometric errors less than 0.04 in the C field: 
(a) in the color-color diagram, and (b) in the color-magnitude diagram. The filled squares are for red giants and asymptotic giants spectroscopically observed by 
Mateo et al. (1998). The filled triangles are for carbon stars spectroscopically confirmed by Azzopardi, Lequeux, & Westerlund (1985, 1986), and the inverted 
filled triangles are for carbon stars photometrically identified by Demers & Battinelli (2002). 
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FIG. 7. — Distribution on the sky of stars selected as "Leo I-like" giants by our color-color and color-magnitude selection methods. The lower panel is a 
blow-up of the central part in the upper panel. The region enclosed by the solid lines is our survey area. The ellipse shows the tidal boundary derived by IH95 
and the hatched region in the bottom panel is outside of our survey area. 
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FIG. 8. — Sky distribution of Leo I giant candidates with the new King limiting radius (solid line) derived in this study. The inner ellipse drawn with a dashed 
line corresponds to the point where the radial density profile starts to deviate from a single-component King model (§4.1). The hatched region is outside of our 
survey area. 
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FIG. 9. — Fits of the stellar distribution of Leo I using (a) a King model with parameters derived by IH95, (b) our new King model with one power-law, (c) new 
King model with two power-laws, and (d) a power-law with core (PLC) model. The horizontal dashed lines are the derived levels of the backgrounds, which were 
subtracted from the data and the model curves shown. The data points in each panel have slightly different positions because the ellipticities and position angles 
derived from the fits require different binning of the data points. The dotted lines in the two mid panels show power laws of the form Sjv oc r~"< , where 7 = 4.4 
for (b), and 7 = 1.7 (shallower slope) and 10.5 (steeper slope) for (c). The errors for each point includes both la Poissonian and background estimation errors. 
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Fig. 10.— 



Radial surface density profiles for the eastern and western halves of Leo I. The closed and open circles are for eastern and western bins, respectively. 



34 



o 

Si 

(0 

o 
Q 

< 



20 



10 



-10 



-20 



— i — i — i — | — i — i — i — | — i — r 

(a) 




i /i /i / 1 . 




20 10 -10 -20 20 10 -10 -20 

ARA (arcmin) ARA (arcmin) 



FIG. 1 1. — (a) Isodensity contour plots of Leo I constructed using Leo I giant candidates. This figure shows the 43 X 43 arcmin 2 region centered on Leo I. 
The newly derived center is at (ARA, ADec) = (0,0). The hatched region is outside of our survey region. The contour levels are 1, 2, 5, 10, 20, 35, 60, 90, 120, 
and 160 stars per pixel, (b) Smoothed isodensity contour plot of Leo I using the same technique for making (a) but with each star represented by a 6f0 wide 
two-dimensional Gaussian (see §4.3). Each pixel has a dimension of roughly O'.lx O'.l. The contour levels are 0.8, 2, 5, 10, 20, 35, 60, 90, 130, and 160 per 
pixel. 
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FIG. 12. — Spectra of four stars observed with the Keck DEIMOS spectrograph. The wavelength region 8450-8750A shown highlights the Ca triplet lines. 
The top spectrum is that of our chosen radial velocity template star (see text), and the bottom three stars are Leo I giants of different brightnesses. 
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FIG. 13. — Telluric offsets for each stellar RV measurement as a function of the position in the mask for (a) KD1 and (b) KD2. The variation in astrometric 
solution from chip to chip in the original Mosaic camera data is revealed by the trends in telluric offset by chip number. 
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FIG. 14. — Histogram for heliocentric radial velocities of 125 stars observed with Keck DEIMOS that have Q>4. Bin size is 5 km s 
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(M - T a ) 

FIG. 15. — Positions of the stars observed with the Keck DEIMOS: (a) in the color-color diagram; (b) color-magnitude diagram. The selection boxes for the 
Leo I giant candidates (from Figs. 4 and 5) are shown in both plots. The filled circles are for stars that have been confirmed to be members of Leo I based on 
their heliocentric radial velocities, while the crosses are for stars confirmed not be members. 



39 




10.16 10.15 10.14 10.13 



RA [h, J2000] 



FIG. 16. — Stars observed with the Keck DEIMOS plotted over the sky distribution of Leo I giant candidates. The approximate field of view for the Keck 
DEIMOS fields are shown with long dashed lines. The inner ellipse shown in solid line represent the King core radius, and the outer two ellipses are same as 
those in Figure 8. The symbols are same as those in Figure 15. 
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FIG. 17. — Distributions of (a) individual heliocentric radial velocities and (b) mean radial velocities along the projected major-axis radial distances. Filled 
squares in the upper panel are stars observed with the Keck DEIMOS, while filled triangles are those from Mateo et al. (1998). 
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FIG. 18. — Radial variation of individual radial velocities and velocity dispersions using (a) elliptical and (b) circular radii. The number of stars for calculating 
each data point (<r v ) in the mid panels are 14, except for the last data point where 11 stars were used. Similarly, 10 stars were used for calculating each data 
point in the lower panels, except for the last data point where 7 stars were used. The vertical dashed lines denote the King limiting radius (geometric mean, 
i.e., nim,geom = HimVl—e for the right panels). The radial velocity outliers (see text) are marked with large circles in the upper panels, and velocity dispersions 
computed without the outliers are shown in open circles in the mid and lower panels. 
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FIG. 19. — Histogram for heliocentric radial velocities of (a) Leo I giants with r ma j/r Um < 0.6, (b) Leo I giants with r ma jjri im > 0.6, and (c) all Leo I giants. 
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FIG. 20. — Same as Figure 16, but the radial velocity outliers in Figure 18 marked in open circles. The ellipse drawn in dotted line corresponds to r ma j/ri im = 0.6. 
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FIG. 21. — Normalized number densities of "extratidal" stars as a function of azimuthal sectors (18° wide) around Leo I. The histogram represents the number 
density in each elliptical sector outside the break radius (major axis radial distance of 10.2 arcmin) — see Figure 8. The final results were divided by the maximum 
number density to scale the relative number densities from to 1. The repeated dashed line at PA = 45° is added for guidance. The broad peaks to the east and 
west reflect apparent tidal arms stretching out from the main body. The excess of density at 180-120° is due to the "bridge" of stars that extends to the southeast, 
and is likely a statistical anomaly. 
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FIG. 22. — Comparisons of observed (left panels) surface radial density profile (upper panels) and velocity dispersion profile (lower panels) with those of 
W-body models. The closed triangles in the model surface density profiles are for the entire sample, while open triangles are for the bound stars. The open circles 
and open squares in the top panels are for unbound stars from the first and second mass loss events, respectively. The symbols used for the velocity dispersion 
profiles are same as those in the mid left panel of Figure 1 8, while the solid and broken lines show the model velocity dispersion profiles for the entire sample 
and the bound stars, respectively. 
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FIG. 23. — Orbital properties of model 117. This orbit is very similar to the one for model 111 since the primary difference between them is the closer 
perigalacticon distance in model 117. White dots (colored Yellow in electronic edition) represent particles that are still bound. Dark gray (magenta in electronic 
edition) and light gray (cyan in electronic edition) mark particles that became unbound in the first and second perigalacticon passages, respectively. The current 
position of M31 is indicated as arrows. 
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FIG. 24. — Comparisons of observed (left panel) surface radial density profile (upper panels) and velocity dispersion profile (lower panels) with those of model 
122. The symbols and lines are same as those in Figure 22 
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FIG. 25 . — The observed Galactocentric radial velocity distribution (leftmost panel) compared with the model velocity distributions using the DEIMOS footprint 
for Leo I. In each model panel, the apogalacticon and perigalacticon distances are noted. The rightmost panel is for a model with opposite orbital pole. 
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FIG. 26. — Comparisons of observed {left panel) surface density profile and velocity dispersion profile with those of models having the same orbital and 
structural properties, but where one has had 2 perigalacticon passage (mid panels) and the other only one (right panels). The symbols and lines are same as those 
in Figure 22. 



